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The tools developed in a preceding article for interpreting spacetime geometry in terms 
of all possible space-plus-time splitting approaches are applied to circular orbits in some 
familiar stationary axisymmetric spacetimes. This helps give a more intuitive picture 
of their rotational features including spin precession effects, and puts related work of 
Abramowicz, de Felice, and others on circular orbits in black hole spacetimes into a 
more general context. 
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1. Introduction 

"Rotating spacetimes" have captured people's imaginations ever since "rigid" 
rotations in Minkowski spacetime were considered within the theory of general 
relativity. Even this simple example which is the foundation of the "fictitious" 
centrifugal and Coriolis forces in classical physics has led to its share of confusion 
about rotation in relativity. Godel's discoverytl of the spacetime which bears his 
name certainly added fuel to the fire, which was again stoked by the discovery of 
the rotating black hole solution of Kerm and its generalizations.tlQ 



The language of gravitoelectromagnetism,El specialized in the preceding com- 
panion articles (to be referred to here as [BCJl]) for stationary axially symmetric 
spacetimes, helps us to understand the effects of rotation as well as those of accel- 
eration and spatial curvature in these three classic spacetime examples. Indeed the 
lines of force of the various gravitoelectromagnetic vector fields, especially in the 
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Kerr spacetimes, help give a more tangible way of interpreting the behavior of test 
particle motions in the gravitational field of these spacetimes. 

Here we focus on the simpler case of circular orbits following Killing trajecto- 
ries in these^ spacetimes, confining our attention to the equatorial plane in the Kerr 
spacetimes .13' lI This test particle motion is an example of "purely transverse" rela- 
tive acceleration.0 By exploring the roles played by the radial spatial gravitational 
forces, one obtains a clearer picture of the action and interrelationships of the var- 
ious gravitoelectric (GE), gravitomagnetic (GM), and space curvature (SC) forces 
that one may define within each spacetime as well as of the correspondences one 
may establish between these different spacetimes. Presenting the details of these 
applications also helps make more concrete the somewhat abstract but powerful 
language of gravitoelectromagnetism itself, which can be valuable in interpreting 
the geometry of other spacetimes. In particular, previous discussions_D|Lcircular 
orbits in black hole spacetimes by Abramowicz, de Felice, and otherstlrEil are fit 
into a more general picture which helps to clarify their particular analyses of the 
behavior of certain properties of these orbits. 

Each of these three classes of spacetimes have natural stationary axisymmet- 
ric nonlinear reference frames, i.e., threaded slicings (hypersurface foliations with 
transversal congruences of curves) of the spacetime which are adapted to two Killing 
vector fields associated with a 2-dimensional stationary axisymmetry group. The 
nonlinear reference frames for the rotating Minkowski and Godel spacetimes have 
an additional translational symmetry making them cylindrically symmetric as well. 
In each case the threading and slicing families of test observers associated with 
these nonlinear reference frames are tied to the geometry of the spacetime and help 
elucidate its properties. 

In the case of Kerr, the nonlinear reference frame associated with Boyer-Lind- 
quist coordinates {t, r, 0, has as its threading observers (following the time 
coordinate line Killing trajectories) the distantly nonrotating observers or static 
observers, while the slicing observers (moving normal to the time coordinate hy- 
persurfaces) are the locally nonrotating observers or zero-angular-momentum ob- 
servers. Both observer families are accelerated, the threading observers opposing 
the dragging along action of the rotating black hole, while the slicing observers are 
dragged along by the hole with respect to spatial infinity. In the Godel spacetime, 
the nonlinear reference frame of cylindrical coordinates {t, p, 0, z} has the thread- 
ing observers moving along the geodesic flow lines of the rotating dust source, while 
the accelerated slicing observers oppose the global rotation of the spacetime. In the 
rotating Minkowski spacetime, the accelerated threading observers are uniformly ro- 
tating, while the geodesic slicing observers are a global family of inertial observers 
associated with the usual time lines of an orthonormal Cartesian coordinate system. 
Each of these three spacetime examples exhibits different configurations of the vari- 
ous gravitoelectromagnetic fields whose comparison offers insights about the nature 
of the spacetimes themselves. 

The Godel spacetime was first studied as an example of a constant gravitomag- 
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netic field configuration by Wilkins and Jacobs,Ea_bascd on the analogy between 
linearized general relativity and electroniagnetisni.L3 Their pioneering work stim- 
ulated the present program of trying to better understand fully nonlinear general 
relativity in terms of observer-splittings. 

2. Gravitoelectromagnetic potentials, fields, and forces 

The "rotating Minkowski" spacetime is just Minkowski spacetime expressed in 
terms of the nonlinear reference frame associated with a family of uniformly rotating 
test observers (threading) and the time hypersurfaces (slicing) associated with the 
global inertial Cartesian coordinates with respect to which this rotation takes place. 
It is most easily described in terms of rotating cylindrical coordinates {t, p,</), z}, 
where 4> = (j)' — Qt gives the relationship to the nonrotating angular coordinate 



angular velocity fl about the z-axis may take any real value (define TZ = l/|f2|), here 
it will be assumed to be positive in order to discuss corotation and counter-rotation 
with respect to the sense defined by the positive 2;-axis and the righthand rule. 

The Godel spacetime is a solution of the Einstein equations with constant dust 
energy density p(Q) and cosmological constant A related hy A — —il,'^ — —(27?.^)^^ — 
— 47rp(o), where TZ = \\/2fl\^^ is Godel's curvature parameter "a" and fl is the 
nonzero constant parametrizing the constant vorticity of the fluid source. Although 
it may take any nonzero value, it will be assumed to be positive here for the same 
reason as above, so that the threading observers corotate with respect to the pos- 
itive z-axis in the cylindrical coordinates used here. These coordinates differ from 
others in the literature by a rescaling, chosen so that to linear order in the vortic- 
ity parameter the Godel and rotating Minkowski metrics are the same, although 
important differences occur at second order. In these coordinates the limiting be- 
havior approaching the axis of symmetry of the associated nonlinear reference frame 
is also the same to lowest order in this parameter, and the limit of both metrics as 
^ is just Minkowski spacetime expressed in nonrotating cylindrical coordinates. 
Thus the new parametrization of the Godel metric emphasizes the key feature of 
its geometry, namely the global rotation, rather than the related spatial curvature. 

The Kerr spacetime describes the geometry around a rotating black hole with 
mass M (directly defining a length scale TZ) and angular momentum per unit mass 
a as seen from infinity. The angular velocity parameter a will be assumed positive 
as above, so that all three spacetime examples rotate in the same sense. The limit 
a = gives the nonrotating Schwarzschild spacetime, while a = is the "extreme 
Kerr" case. The usual Boyer-Lindquist coordinates {t, r, 6, (j)} will be used here. 

In each spacetime the Killing vector field cq" = S^t generates the stationary 
symmetry while the spacelike Killing vector field S"^ generates the axisymmetry, 
and the spatial coordinates are orthogonal. For uniformity of discussion, the radial 
Kerr coordinate symbol r will be used to denote the cylindrical coordinate p in 
the other two examples as well, while the physical (orthonormal) component along 
—S^g perpendicular to the equatorial plane in Kerr will be referred to as along the 




Although the coordinate 
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positive z-axis as in the case of the plane z = in the other spacetimcs and will 
be indicated by the index z. The latter plane will be referred to as the equatorial 
plane in all cases. For graphing purposes it is natural to introduce the rescaled 
radial variable f = r/TZ for all these spacetimes and a = a/TZ = a/M for the 
additional parameter in the Kerr case. 

For all the three spacetimes under consideration, the spacetime metric restricted 
to the world sheet of this plane can be expressed in the reference (i.e., coordinate), 
threading, and slicing decompositions respectively as 

ds^ = ^'^'^gttdf + 2^^^gt4,dtd(t> + '^'^^g4,4,d(j>^ + ^^''grrdr'^ 
= ~M^{dt- M4,d(l)f +-f^^d(t>^ +-frrdr^ 

= -N^dt^ + g4,4d(j) + N'^dtf + grrdr'^ . (2.1) 

Since the threading and slicing observer-adapted frames associated with this adapt- 
ed coordinate system are obtained by projection of the coordinate frame, their 
spatial structure functions vanish. Table |l| gives the expressions of the various 
threading and slicing quantities for the three cases. 

Table 1. The lapse, and the nonzero observer-adapted components of the shift and spatial metric 
in both the threading and slicing points of view are given for the rotating Minkowski, Godel and 
Kerr spacetimes on the equatorial plane. Note the Godel parameter relation \/2|f^| = 1/7?. useful 
for reinterpreting quantities in terms of spatial curvature rather than vorticity. 



lapse, shift and 
spatial metric 


Rotating 
Minkowski 


Godel 


Kerr 


Ar= (-(4)g00)-l/2 


1 


C/Vl-S2 




N'l' = N^/g^^ 


a 


-s2) 


-2aM/{r^ + a?r + 2o?M) 








-2aM/r 


g<t>4> = '*^9</.</) 


r2 


2s2(l - s2)/n2 


(r^ _|_ a'^r + 2a?M)/r 


(4) 

Qrr — 'Jrr — grr 


1 


1 


r2/A 






1 


^(r - 2M)/r 






2s2/r2 


-2aM/(r - 2M) 


= ("^0 - {'■%o^f/'-^^goo 


7^2 


2s2c2/n2 


rA/{r - 2M) 



7 = (1 - n2r2)-i/2 ^ s = sinh(y2nr/2) , c = cosh(V2Qr/2) , t = tanh(V2nr/2) , 
S = sinh(v^nr) , C = cosh(V2nr) , T = tanh(V2f7r) , A = - 2Mr + , 



The general formulas of [BCJl] for constant speed test particle circular orbits 
in stationary axisymmetric spacetimes are easily evaluated for the present explicit 
metrics. The test particle moves along the direction with constant speed. The 
4-velocity of a nonzero rest mass test particle is parametrized by the coordinate 
angular velocity C = = dc^jdt as follows 

[/" = r[<5"t -f C^"^] , (2.2) 
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where T = dt/ dru > is defined by 

r-^ = + 2C(%^ + c^(%^] = -(%^(C-C-)(C-C+) 

= M^l - M^Cf ~ ^^^C^ = - g^^iC + (2.3) 

and Tij is a proper time parametrization of the world line. The timelike condition 
for the 4-velocity requires F^^ > 0, constraining C, to belong to the interval 
[C,-, (^+] between the roots of the quadratic equation T^^ = in C corresponding to 
null directions, namely 

The 4-velocity of a zero rest mass particle (for which F"^ = 0) has an arbitrary 
normalization factor r(,iun) in place of F in equation (2.2) 



P±=r^nun)[S''t + (:±S%] . (2.5) 

Circular orbits for which C > or C < will be referred to respectively as corotating 
or counter-rotating (with respect to the nonlinear reference frame or the threading 
observers) . 

Note that the coordinate angular velocity of the slicing observers is just the 
average of the two limiting angular velocities 

C(sl) = C(nmp) = (C- + C+)/2 = -N't" , (2.6) 

which is just equation (33.16) of Misner, Thorne, and Wheeleiil in the specific 
context £if the Kerr spacetime. Their exercise (33.3) following the discussion of 
Bardeenc3 gmplies to the general (orthogonally transitive) stationary axially sym- 
metric case,c3 so that one may interpret the slicing observers as the locally nonro- 
tating observers with respect to the Sagnac effect. They experience no Sagnac effect 
for the oppositely directed accelerated photons constrained by mirrors or fiber opti- 
cal cable to remain on a given circular orbit, meaning that the alternating meeting 
points of these photons lie on the same observer world line. They are also called 
the "zero angular momentum observers" (ZAMO's) since they are orthogonal to the 
angular Killing vector and therefore have vanishing angular momentum. A comple- 
mentary formula exists for the angular coordinate component of the threading shift 
1-form 

M^ = {(:--'+C+-')/2. (2.7) 

This is related to the Sagnac effect as well, as explained below. 

The physical components of the velocities measured by the threading and slicing 
observers for such motion are related to the coordinate angular velocity by linear 
or fractional linear transformations 

iy{U,m)^ = 70//'C/[M(l - A/^C)] , J^(C/, n)^ = g^//2(C + 7V0)/iV (2.8) 
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and 

C = Mv{U,mf/[-i^^''''' + MM^v{U,mf] 

= -N't" + Ng^^-^'^u{U,nf . (2.9) 

Note that when the shift is nonzero, test particle motions with angular velocities of 
equal magnitude but opposite sign lead to physical velocities which do not have the 
same magnitude and vice versa. When i'{U,u)'^ = ±1, the latter equation reduces 



to Eq. (2.4) 



The "coordinate" gamma factor is easily expressed in terms of the usual Lorentz 
gamma factor associated with these relative velocities 

r = 7(C/,m)/[M(l-M^C)]=r(t/,m) 

= "f{U,n)/N = T{U,n) . (2.10) 

These formulas may be used to express the angular momentum (per unit mass) 

P^^U^^ Vv{U,n)^ = ff00r(C- C(sii)) (2.11) 

of U defined by the rotational Killing vector J"^ and its Killing energy (per unit 
mass) £ = —Ut = M~^j{U, m) defined by the Killing vector S^t, both conserved for 
geodesic motion. These are related to the coordinate gamma factor by the identity 
-1 = = T{-£ + (p^) in the timelike case and = P„P" = r(„uii) + CP4,) 

in the null case, where £ = —Pt- Note that the slicing relative velocity is directly 
proportional to the angular momentum. 
In the timelike case the ratio 

^^^^_5« + Cff^ (2.12) 

P<l> 9t<p + (,94,4, 

defines the angular velocity of the spacelike circular orbit orthogonal to [/" with 
unit tangent ~ r((5"t + C'^"^), T > and having the same sense of rotation. 
This is the angular direction of the local rest space of the test particle. 

Consider only a nonzero rest mass test particle in what follows. The various 
spatial forces acting on such a particle all point along the radial direction whether 
expressed in the threading or in the slicing observer-adapted frame. In the "spatial 
equation of motion" (9.9) of [BCJl], the Lie total spatial covariant derivative of the 
spatial momentum reduces to minus the space curvature force (see equations (12.7) 
and (12.23) of [BCJl]) and one finds the following simple results for the equation 
of motion in the threading, hypersurface, and slicing points of view respectively 

-F{U,uf = Fg(C/,i.)^ + F(sc)([/,„)^ , u = m,n, 
-F{U,nY = F(g(C/,n,eo)^ + F(s^)(C/,n,eo)' , (2.13) 

namely, minus the relative non-gravitational spatial force must balance the sum of 
the Lie spatial gravitational force and the space curvature force. The spatial grav- 
itational forces in the various points of view can be separated into their gravito- 
electric (GE), vector gravitomagnetic (CM), and symmetric tensor gravitomagnetic 
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expansion (EX) components 

Fg}(C/, my = -f{U, m)[g{my + v{U, mfHimy] 

= F^-^'^^U^my + F(-'^^\U,my , 

Fl^^{U,ny = j{U,n)[g{ny-2,.{U,nye{ny^] (2.14) 

= F(G'^)(C/,n)^ + F(^^)(C/,r^)^ 

43(t/,n,eo)^ = ^iU,nMny + \v{U,nyH{n,eoy - y{U,nfe{ny^] 
= ([/, ny + F(g^) ([/, n, eo)^ + F^^^^) ([/, n, cq)^ , 

where the sum of terms defines respectively the individually named force terms. 
The space curvature forces reduce simply to the sign-reversal of '){U, u) times the 
Lie relative centripetal acceleration 

F^^'^\U,my = -K{(t>,my-f{U,m)\i^{U,my\'^ , 

F^^^\U,ny = -K{(j>,nyj{U,n)\u{U,n)^f , (2.15) 
F^^^\U,n,eoy = -«(</., n)^7(f/, ri)K[/, n)^]K«7, n)^ - Keo,n)"^] , 

where 

K{(P,my = -(In 7^^1/2) . ^ K(0,n)^ = -(In 3^^1/2),^ (2.16) 

are the signed Lie relative curvatures of the ^ coordinate lines in the threading and 
hypersurface points of view and 

u{eo,n)^ = N-^N^ = {g^^f'^N-^N'^ (2.17) 

is the physical component u{m,ny along the positive (j) direction of the relative 
velocity of the threading observers with respect to the slicing observers when eo" 
is timelike, and of the relative velocity of the time coordinate lines in general. 
The caret index notation indicates physical (orthonormal) components along the 
orthogonal coordinate frame vectors in the slicing point of view and along their 
corresponding spatially projected vectors of the associated orthogonal observer- 
adapted spatial frame in the threading point of view. 

Finally, the nonzero components of the gravitoelectric and gravitomagnetic vec- 
tor fields together with the expansion can be calculated from their expressions in 
the observer adapted frames (see equations (14.5) and (14.6) of [BCJl]) 

g(m)r = {lTrY''^9{m)f = -(InM),^ , 

9{n)r = {grry^'^g{n)r = -{lnN)^r , 

H{my = H{my = M{^„l^^)-^'^M^^r , 

H{n,eoy = H{n,eoy ^ N-\grrg4>4>)~^^^N^,r , 

0{ny^ = {9^^9rrY'^0{n),^ = -{l/2)g^^N-^N\r ■ (2.18) 
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Table 2. Causal restrictions, spatial gravitational force field quantities, and geodesic and null 
conditions for circular orbits in the rotating Minkowski and Godel spacetimes. 





Rotating 
Minkowski 


Godel 


length scale 7^ 




1 V " 1 


threading region 
of validity 


r < r(h) = 


everywhere 


slicing region 
of validity 


everywhere 


r < r(ij) = 7?.arcsinhl 


u(m,n)'^ = —u(n,m)'^ 


fir 


V2t 


gravitoelectric 
field 


g{mY = -y'-'n'^r > 
9{nf = 


g(mY = 


gravitomagnetic 


H{m)^ = 2^7^ > 


H{mY = 2Q > 


field 


H{n,eoY = 2n 


H{n,eoY = 2C> 


expansion 


Hn),^ = 


e(n).7 = < 


signed relative 


K{(f,,mY = < 




curvature 


K.{4>,n)f- = - i < 




counter-rotating 


</>_ = -n 




and co-rotating 
timelike geodesies 


u(U-,mY = - fir 
u(U-,nY = 


v{U±,mY = - V^T ,0 
v{U±,nY = -f^ ,y2t 


null orbits 


C± = -n± 1/r 


C± = tS["1±*"VV2] 



Explicit expressions for all of these quantities in the various points of view are given 
in Tables || and || for each of the three spacetimes under consideration. 

The threading point of view is valid for those points of spacetime where the 
threading is timelike, i.e., where '-^■'gtt < (or M > 0), while the hypersur- 
face/slicing points of view are valid where the slicing is spacelike, or equivalently 
where the normal to the slicing is timelike, i.e., where '•''^g** < (or > 0). If 
the threading observers become spacelike in a region where the slicing observers are 
still timelike, then the threading point of view is still valid where the magnitude of 
the relative velocity of the threading observers with respect to the slicing observers 
is \v{m,n)'^\ = |A^~^iV'^| < 1. Vice versa, if the slicing observers become spacelike 
in a region where the threading observers are still timelike, then the slicing point of 
view is still valid in the region where \v{n^ m)'^\ < 1. The regions of validity of the 
threading and slicing points of view are given in Tables | and |. 

For both the rotating Minkowski and Godel spacetimes, the magnitude of the 
relative velocity is an increasing function of r which leads to the existence of an 
outer light cylinder at r = r^h) ("h" for observer "horizon") where the worldlines 
of one of the families of test observers become null and the corresponding point of 
view is no longer valid. In the rotating Minkowski case the light cylinder exists only 
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Table 3. Causal restrictions, spatial gravitational force field quantities, and geodesic and null 
conditions for circular orbits in the Kerr spacetime. 





Kerr 


length scale 


A4 


threading region 

of validity 


■>■ > '■(erg) = 2>1 


slicing region 
of vuliHitv 




u{m,n)^ = —u{n,m)^ 


2aM 


gravitoelectric 
field 


n(mY = < 


gravitomagnetic 




field 


H{n,eor = ^>0 


expansion 




signed relative 
curvature 




counter-rotating 


0+ 7—^ 

l±ay'M/r'> 


and co-rotating 


u(U±,m)^ = ^^^^-^ 

^ ' a±(r-2M)^r/A4 


timelike geodesies 


uiU±,n)i' - o^y^'^^+r^ 


null orbits 





in the threading point of view and it occurs at the radius for which \v{m, n)'^\ = 1, 
while in the Godel case it exists only in the slicing point of view where it occurs 
when \i'{n,m)'^\ = 1. However, beyond that horizon g^^ becomes negative, leading 
to the famous closed timelike (j) coordinate lines. In the Kerr case where discussion is 
confined to the equatorial plane, the situation is reversed and the relative velocity is 
a decreasing function of r leading to an inner light radius. The slicing point of view 
is valid outside the event horizon which occurs at the value of r for which A = 0, 
while the threading point of view is valid for r > r"(ei.g) outside the ergosphere (which 
in turn surrounds the event horizon), where Tf^^gj is determined by the condition 
n)'^\ = 1. The Godel slicing observers attempt to resist the global rotation of 
the spacetime, but are forced to corotate at the outer observer horizon; similarly 
the threading and slicing observers in Kerr are forced to corotate at their respective 
inner observer horizons, namely the ergosphere and the event horizon. 

The gravitoelectric field is the sign-reversed acceleration of the observer congru- 
ence, thus revealing the accelerations which characterize the threading and slicing 
observers. For example, in the rotating Minkowski spacetime the threading ob- 
servers are accelerated radially inward while the counter-rotating slicing observers 
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(at rest in a global inertial frame) have zero acceleration. In the Godel spacetime 
the situation is reversed and the threading observers arc not accelerated while the 
counter-rotating slicing observers are accelerated radially outward. In both cases 
an outward acceleration must be added to resist the global rotation of spacetime 
(more precisely, of the nonlinear reference frame) by counter-rotating. The equato- 
rial Kerr slicing and threading observers are both accelerated radially outward to 
oppose the attraction of the central mass, leading to inward gravitoelectric fields 
which allow circular orbits even in the nonrotating Schwarzschild limit. In the Kerr 
case the threading observers counter-rotate with respect to the slicing observers 
and one finds that their acceleration is larger then the acceleration of the slicing 
observers. 

In each case considered here the gravitomagnetic fields are along the positive 
z-direction. In both the rotating Minkowski and Godel spacetime this is due to 
the fact that the shift 1-forms are along the positive (j) direction and their physical 
components along that direction are increasing functions of r, while in Kerr space- 
time the same 1-form physical component is a decreasing function of r but its sign 
is reversed. Thus in each case the sign of the radial gravitomagnetic force compo- 
nent depends only on the sign of the relative velocity of the test particle along the 
angular direction, namely a positive (outward) force for co-rotating orbits and a 
negative (inward) force for counter-rotating orbits, where here the terms corotat- 
ing and counter-rotating are with respect to the given observer family. Note also 
that in the Godel spacetime the threading and slicing gravitomagnetic fields are the 
same and uniform (spatially covariant constant). That they are equal can also be 
seen directly from the transformation law for the gravitomagnetic vector field given 
in equation (11.6) of Ref. || and using the fact that the threading gravitoelectric 
field and the Lie derivative along cq" of the shift are zero and that the relative 
projection P{m,n)~^ reduces to the identity along the radial direction orthogonal 
to the plane of the relative motion of the two observers. As discussed in [BCJl], the 
expansion tensor, zero in the threading point of view, has one possibly nonvanishing 
r-cf) component which is zero (rotating Minkowski) or negative (Godel and Kerr) in 
the hypersurface and slicing points of view. Like the radial gravitomagnetic vector 
force, the radial expansion force is also positive (outward) for co-rotating orbits 
and negative (inward) for counter-rotating orbits, again with the sense of rotation 
referred here to the observer family. The sum of the gravitomagnetic vector force 
and the expansion force equals the total gravitomagnetic tensor force as discussed 
in [BCJl]. 

On the other hand, the signed relative curvatures can change sign within the 
range of validity of one point of view yielding radial centripetal accelerations (sign- 
reversal of the space curvature forces) that can be either inward (k(</), uY < 0, the 
usual case) or outward («(</', uY > 0). Figure |l| shows k(0, uY for the various cases 
in the threading and hypersurface points of view. The local extrema of g^^ and 
7^^ are the points where the relative Lie centripetal acceleration changes sign and 
the relative curvature of the spatial trajectory vanishes, yielding the Lie relatively 
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straight trajectories which are also spatial geodesies in each point of view. 

In Figure ^and successive figures, the Godel and Kerr diagrams are divided into 
regions A, B, and C by thick vertical dashed lines according to whether 2, 1, or of 
the two oppositely directed circular geodesies are timelike. The boundaries of these 
regions, as discussed in detail in section 4, occur at the radii of the two null circular 
geodesies. Thin vertical dashed lines also mark the horizon and outer ergosphere 
radii in the Kerr case. 

The behavior of the Lie relative curvature depends on the intrinsic geometry of 
the Riemannian metric on the quotient space appropriate to each point of view. In 
the threading case this metric is just the natural projection of the spatial metric 
to the observer-quotient space, while in the slicing point of view, it is instead the 
pullback of the spatial metric to the quotient by the threading congruence. For the 
case of planar orbits, 2-dimensional embedding diagrams in either 3-dimensional 
Euclidean space or 3-dimensional Minkowski spacetime M3 as appropriate are 
useful to interpret the effects of this spatial geometry on the orbits in this plane, as 
well as on the precession of the spin of a gyroscope following such an orbit. ThorneEJ 
has given an neat illustration of this latter effect in terms of the conical defect of the 
tangent cone in the embedding space. The details of the embedding are explained 
in the appendix, together with the actual diagrams of the cross-sectional curves of 
these surfaces of revolution. 

3. Circular geodesies 

The timelike circular geodesies are the circular orbits along which the rela- 
tive spatial force F{U, uY vanishes. The relative velocities corresponding to these 
geodesies can then be found by setting to zero the sum of the Lie spatial gravita- 
tional force and the space curvature force, leading to a quadratic equation in those 
velocities or in the coordinate angular velocities (linear when K{(f),uy ~ 0). The 
regions where the geodesies are timelike are those regions where the corresponding 
spatial relative speeds are smaller than 1 (provided that the observers are also time- 
like). Let U±., i^{U±,u)°', and 0± be the 4-velocity, relative velocity, and coordinate 
angular velocity of the corotating (+) and counter-rotating (— ) geodesies in each 
spacetime, when they exist. 

For example, in the rotating Minkowski case there are only counter-rotating 
circular geodesies, which correspond to the points fixed in the global inertial frame 
with respect to which the rotating nonlinear reference frame rotates, namely the 
orbits of the slicing observers {v{U-, n)'^ = in Tables |2| and ^. These worldlines 
are clearly timelike everywhere. Their spatial velocity relative to the threading 
observers becomes larger than the velocity of light when ilr > 1 only because the 
threading observer congruence becomes spacelike there. In the Godel spacetime the 
co-rotating geodesies are instead the orbits of the threading observers {v{U+, m)'^ = 
0). These worldlines are then timelike in the region of validity of the threading 
point of view, i.e., everywhere. On the other hand the counter-rotating geodesies 
are timelike as long as \v{U-, m)'^\ < 1 or |j/(/7_, n)''^| < 1 which gives r < ^'''(h)- In 
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Fig. 1. K versus f: plots of the signed Lie relative curvature k(4>,uY of the circular orbits in 
the equatorial plane of the rotating Minkowski, Godel, and a/JVl = 1/2 Kerr spacetimes in the 
threading (u = m) and hypersurface {u = n) points of view as functions of f = r/li. (Minkowski, 
Godel) or f = r/JV[ (Kerr). In this particular Kerr case the relatively straight circle occurs just 
inside the region C (which begins at f 2.35) described below. 
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the Kerr spacetime the equations \h'{Uzp,n)'^\ = 1 and \iy{Uzp,m)'^\ — 1 have roots 
(vahd only within the region of vahdity of each point of view) at r(r — 3M) = 
±2aAi y/rJjA . In the extreme Kerr case a = A4, these roots are r = A4 for the 
co-rotating geodesies and r = AM for the counter-rotating geodesies. 

Tables || and ^ show that for the Godel and Kerr spacetimes which have two 
oppositely directed circular geodesies, both the magnitude of the coordinate an- 
gular velocity C and of the angular component of the velocity i'{U,u)'^ are larger 
for geodesies in the counter-rotating direction than in the corotating direction in 
both the threading and slicing points of view. This immediately implies that the 
geodesic gamma factor j{U±, u) and with some additional reasoning the coordinate 
gamma factor T{U±^u) are also both larger in the counter-rotating direction than 
in the corotating direction. The latter follows from the former using the slicing 



representation of the coordinate gamma factor given in Eq. ( 2.1C ) . This general 
counter-rotation effect manifested in this asymmetry between the co-rotating and 
counter-rotating directions is the direct consequence of an upward gravitomagnetic 
field which adds an inward radial force to the total sum for the counter-rotating case 
increasing the speed, and an outward radial force to the total sum for the corotating 
case decreasing the speed. 

The asymmetry in the coordinate angular velocity between the corotating and 
the counter-rotating circular geodesies is the origin of the precession of their alter- 
nating meeting points (after a full revolution, not half a revolution) in the counter- 
rotating direction. Analogous to the zero Sagnac effect (slicing) observers which 
follow Killing trajectories containing the meeting points of oppositely directed cir- 
cular null paths, one caii^introduce "geodesic meeting point observers" (only in the 
equatorial plane in Kerrcl) containing the meeting points of the oppositely directed 
circular geodesies. Their angular velocity is analogously the average (see Tables | 
and I) 

C(gmp) =(</>-+ </'+)/2 . (3.1) 

For Kerr this has the negative value C(gmp) = —aA4r^^/{l — a'^Air^^), and for 



Godel it is also negative C(gmp) — 4'-/'^ since 0+ = 0. Using Eqs (p.§|), (p.9|), one 
can show that the slicing velocities are also related by averaging 

HUis.-np),nf - [iyiU-,n)^ + iyiU+,n)'f']/2 , (3.2) 

as is trivially the case for the slicing velocities of the null meeting point observers 
and the oppositely directed null paths. 

The spatial equation of motion evaluated along the circular geodesies with zero 
total spatial force describes how the various radial spatial forces balance. In the 
threading point of view the expansion is zero because of the stationary symmetry so 
only the gravitoelectric, the vector gravitomagnetic, and the space curvature forces 
are present. In the hypersurface point of view the gravitomagnetic field vanishes 
and one is left with only the vector gravitoelectric, expansion, and space curvature 
forces. In the slicing point of view all the forces contribute to the radial force balance 
equation. Figures ^ and ^ show the plots of the various radial forces (divided by 
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the common factor ^{U, u)) evaluated along the circular geodesies as functions of 
the radial coordinate r in the several points of view and for the various cases. To 
better visualize the results, some suggestive 3-dimensional diagrams showing the 
balance of the various forces in space may be found in Carini, Bini and Jantzenjj 
with relative magnitudes of the various forces given in the limit of small rotational 
speeds, i.e., near the axis of symmetry in the rotating Minkowski and Godel cases, 
and far from the black hole in the Kerr case. We now discuss the force balance for 
each of the three spacetimes under consideration. 



3.1. Rotating Minkowski spacetime 



In the rotating Minkowski spacetime the radial force equation reduces respec- 
tively in the threading, hypersurface, and slicing points of view to 



-F{U,nY 
'F{U,nY 



|i/(C/,m)'^|2 



-/iU,m) — {i^{U,m)'t' + nrf , 
r 

r 

F^'^'^^iU^n^eoY + F^''^\U,n,eoY 



(3.3) 



Table ^ shows the explicit expressions of the various spatial forces evaluated along 
the circular geodesies where they balance. 



Table 4. The various spatial radial forces divided by a common factor -f{U, u) evaluated along the 
equatorial counter-rotating circular geodesies arc given as functions of r in the various points of 
view. 



Rotating Minkowski Spacetime 






^-l^(GM)f 


^-l^(EX)f 




Thd_ 











Hyp- 














Sli- 















The geodesies are the orbits of particles at rest in the associated global inertial 
frame with respect to which the rotation takes place. In the threading point of 
view these orbits must counter-rotate in order to compensate for the rotation of 
the threading observers and the outward gravitoelectric force and space curvature 
force along them are equal and add together to balance the inward gravitomagnetic 
force. On the other hand the slicing observers are simply the fixed observers in the 
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(b) Minkowski: Tlid 



0881 T76T 

(c) Godel: Sli - 





um\ ^ 

(f) Godel; Hyp + 



1.763 




1.763 



Fig. 2. F versus f: plots of the space curvature (SC), gravitomagnetic (GM), expansion (EX), and 
gravitoelectric (GE) spatial radial forces (divided by a common 7(C/, u) factor and multiplied by 
TV) along the corotating (+) and counter-rotating (— ) circular geodesies in the rotating Minkowski 
and Godel spacetimes in the threading (thd), hypersurface (hyp) and slicing (sli) points of view, 
plotted versus the rescaled radial coordinate f. The horizontal ranges for the Godel counter- 
rotating/corotating cases are respectively r £ [0,r(ij)/2] and r £ [Oi''(h)]- The thick dashed 
vertical line separates the region A near the origin where both circular geodesies are timelike from 
the region B where only the corotating geodesic is timelike. The observer horizons are marked by 
a thin dashed vertical line. 
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(a) Kerr; Thd - (b) Kerr: Thd + 

0.2 0.2 




Fig. 3. F versus f: plots of the space curvature (SC), gravitomagnetic (GM), expansion (EX) and 
gravitoelectric (GE) spatial radial forces (multiplied by J\4/f{U,u)) along the corotatlng (+) and 
counter-rotating (— ) equatorial circular geodesies around a Kerr black hole with a/Ai = 0.5, In 
the threading (thd), hypersurface (hyp), and slicing (sll) points of view as a function of f = r/Ai. 
The thick dashed vertical lines at f 2.35,3.53 separate the three regions A which extends out 
to infinity, C which extends inward to the horizon (f(hor) ~ 1.87) and B which is sandwiched in 
between them and which contains the ergosphere boundary (f{crg) = 2). The observer horizons 
are marked by thin dashed vertical lines. 
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nonrotating Minkowski spacctime and so they experience no spatial gravitational 
field in the hypersurface point of view. However, in cylindrical coordinates they 
measure a radially outward space curvature force, which vanishes only along the 
circular geodesies and it is due to the curvature of the circular coordinates rather 
than to a real curvature of the space which is actually flat. 

In the slicing point of view there is a gravitomagnetic force and a space curvature 
force but again both vanish along the circular geodesies. The relation between the 
threading and the hypersurface points of view in this case is given simply by the 
addition of velocity formula 



1-^ - 



v{U,mY + Vlr 
-i-f{U,m){l + v{U,m)^nr) , 



(3.4) 



liU,n) 

which just shows that the two points of view are related by a boost. 

3.2. Gddel spacetime 

The radial force equations for the circular orbits in the Godel spacetime in the 
threading, hypersurface, and slicing points of view respectively are 



F{U, m) 



'F{U,nY 



^/2nT-^\v{U,mf\ 



2VLv{U,m)'l' 



2t(l-s2) 



'F{U,nY 



1 - s 

i^(SC)(C/,n,eo)^ + F(«^)(t/,n,eo 
+F(GE'(;7,n)^ 

f [^/2^7(l-2s2) 

7(C/,7i) 



niy{U,n)'t' 



2t(l-s2) 



1 - s 



F(^^Hu,n,eor 



1 -S2 



(3.5) 



Table || shows the values of the various spatial forces evaluated along the two sets 
of geodesies where the forces balance. 

In this case the threading observers move along (corotating) geodesies which 
are the same trajectories as the dust particles, while a second family of geodesies 
counter-rotates as in the rotating Minkowski case. In the threading point of view 
the outward space curvature force balances the inward gravitomagnetic force along 
these counter-rotating geodesies, while all forces vanish for the corotating geodesies. 
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Table 5. The various spatial radial forces divided by a common factor 7(C/, u) evaluated along 
the counter-rotating (subscript — ) and corotating (subscript +) circular geodesies are given as 
functions of r in various points of view. 



Godol spacetimo 




^-l^(SC)f 


^-l^{GM)f 






Thd_ 


2V2m 


-2v/2r2T 








Hyp- 


V2nt 





2^/2(1 ts^ 






{l-s2){l-232) 


1-3^ 


sn_ 












l-2s^ 


(1~S^){1~23^) 


1-3^ 


Thd+ 














Hyp+ 


V2rit(i-23^) 





2v^nts^ 

1-3^ 


1-3^ 


Sli+ 





V2nt 


1-3^ 


1-3^ 



In the slicing point of view, the space curvature force for the corotating geodesies 
vanishes due to the factor [i^{U,n)'^-N-^N'>']i'{U,n)'t' = N-'^^iy{U,n)'^ in equation 
(2.15) which vanishes for (jj = 0. Both the hypersurface and shcing space curvature 
forces for the counter-rotating geodesies and the hypersurface space curvature force 
become negative, corresponding to a Lie centripetal acceleration directed radially 
outward, for r > r(^^s) > ^(h)/2, where r(i.s) = 2arcsinh(2~-'^/^) « 1.317 defines the 
Lie relatively straight trajectory. This effect is closely related to the result discussed 
at length by Abramowicz et al in the case of circular orbits in the Schwarzschild 
spacetime and just corresponds to the turning back of the embedded 2-surface of 
revolution towards its axis as shown in the appendix. However, the counter-rotating 
geodesic is only timelike for r < r(h)/2. 



3.3. Kerr spacetime 

Here one finds the following spatial forces respectively in the threading, hyper- 
surface, and slicing points of view 



F{U,m) 



'r{r-2Mf -Ma^ 



j{U, m) 



+ 



2aM 



r2VA(r - 2M) 



-v{U,mf 



\v{U,m) 



MVA 



r'^{r - 2M) 

i^(SC) (^^ ^y + i.(EX) ^y ^ p(GE) ^ 

J\ \/A(r3-a2M) 



r2(r3 + a^r + 20^]^) 
2aX(3r2 + a2) ^ 



r^{r-2M) 



r^[r^ -f a?r + 2a?M) 
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-FiU,nY 



r^VA{r^ + a'^r + 2a'^M) _ 



2aM(a^M-.3) ^^^^^^^ 



7-3 (r3 + a^r + 2a'^M) 



— „{U,nf 
M[(r2+a2) 



a7W(3r2 



r2(r3 + a^r + 2a'^M) 



v{U, n) 



- Aa^Mr 



(3.6) 



In both the threading and shcing points of view the addition of an upward (+z di- 
rection) gravitomagnetic field, in comparison with the corresponding Schwarzschild 
case, adds a radially inward gravitomagnetic force for the counter-rotating circular 
geodesies and a radially outward one for the corotating geodesies. The coordinate 
angular velocity of the counter-rotating geodesies thus increases in magnitude 
with respect to the Schwarzschild case while the coordinate angular velocity 0+ of 
the corotating geodesies decreases in magnitude. 

In the threading point of view, the space curvature force changes sign at the 
single real root r(j.s) of the equation r(r — 2AAY — Ma? — Q which occurs outside 
the ergosphere in regions C and B, as shown in Figure ^, thus only affecting the 
corotating geodesies in region B. This corresponds to a Lie relatively straight tra- 
jectory which in the extreme case a = 1 occurs at fu^^ = (3 -I- \/5)/2 w 2.618, 



while for the example a = 0.5 one has r. 



(rs) 

(rs) ~ 2.328. Note that the Lie centripetal 
acceleration points radially outward for r < r(j.s) while in Godel it points outward 
for r > r(i.s), i.e., approaching the observer horizon in both cases. 



4. Accelerated circular orbits 

The physical non-gravitational force acting on a test particle that moves with 
4-velocity U is just the 4-force /(J7)" which is related to the rescaled apparent 
3-force by the equation (9.9) of [BCJl], namely 



Fiu,ur^jiu,u)-'p{urf3f{uy 



(4.1) 



In the purely transverse relative acceleration applications that are considered here, 
the 4-force is along the radial direction where the projection reduces to the identity. 
The equation of motion may then be simply expressed in terms of the physical force 
by equation (14.7) of [BCJl], namely 

fiUy = jiU,u)F{U,uf ^ j{U,u)[-F<^''^\U,ur - fI^^{U,uY] = aiUy . (4.2) 
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Kerr: r versus a/m 

4 




0.2 0.4 0.6 0.8 



Fig. 4. f versus a: a plot of the Kerr Boyer-Lindquist radial coordinate f versus d = a/M £ [0, 1] 
for the horizon (hor), outer ergosphere boundary (erg), the two boundaries between the regions 
A, B, and C, and the threading relatively straight circles (rs) and embedding signature-switching 
circles (ss) for the equatorial plane. 



The middle equality is just the first of equations ( ^.13 ). This relationship enables 



one to express the single observer-independent quantity in terms of all the vari- 
ous points of view using equations (2.14) and ( ^.15 ). However, it is important to 



remember that for relative motion not of this kind, the situation is much more 
complicated. 

In the present case the expression for the physical force (namely, the 4-accelera- 
tion by the equation of motion) is a quadratic expression in the relative velocity 
physical component i'{U,u)'^ multiplied by the square of the associated gamma 
factor, leading to a quotient of two quadratic expressions multiplied by a space 
curvature factor. The zeros of the denominator occur when the velocity approaches 
±1, i.e., corresponding to photon orbits, values which bound the domain (—1, 1) of 
allowed values of the velocity. The zeros of the numerator are just the velocities 
which correspond to circular geodesic motion and may occur inside (subluminal) or 
outside (superluminal) the physical domain (— IJ.) of the function of the velocity. 

As noticed by Barrabes, Boisseau and IsraelEJ in their discussion of the hyper- 
surface point of view for Kerr and Schwarzschild, it is very useful to express the 
numerator of the physical radial force expression in the explicitly factorizcd form 

fiUy - Kic^,uYjiU,ufHU,u)'^-^iU.,u)^MU,u)^-,.{U+,u)'^] 



(C-C^)(C-C+) 



(4.3) 



valid when k(0, uY ^ 0, where the roots of these factors in the first case are the 
geodesic velocities {v{U-,uY < i^{U+,uY), at least one of which blows up when 
k(0, uY = so that the force remains finite when other terms besides the centripetal 
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acceleration term are present. This also applies to the Minkowski spacetime with 
v{U-,u)^ = i^(C/+,u)'^. Although the right hand side of this equation contains 
observer-dependent quantities, the expression itself is observer-independent since it 
represents the physical force acting on the test particle. Thus in general, even when 
the spatial curvature force changes sign due to a change in sign of u)*", the total 
force will not. The analogous expression for the zero rest mass particle is 

!{P±Y = K(cl,,uYE{P±,uf[±l-iy{U-,u)'f'][±l-,y{U+,uf] 

- 4i?(P±,n)2«;(0,nf(C±-0_)(C±-(/-+)/(C--C+)' • (4.4) 



de FeliceEZlEEl first used such a factorization (4.3) of the physical radial force 
for Kerr expressed in terms of the coordinate angular velocities rather than the 
physical component of the relative velocities along the orbit. The new angular 
velocity parameter y — ^/(l — a(^) (a fractional linear tranformation) he uses in 



place of the parameter ^ (his il) of equation (2.2) for constant angular velocity 
world lines corresponds locally but not globally to a new choice of slicing by a local 
time function ±!_ = t — acf) as discussed in the appendix of Greene, Schucking, and 
Vishveshwara.c3 Although a global such function does not exist, one can introduce 
it as a well-defined change of parametrization on the world lines of test particles. 
The advantage of this is that the new expressions for the coordinate angular velocity 
d(j)/dt' of timelike circular geodesies in the equatorial plane of Kerr are reduced to 
the simpler Schwarzschild form, while retaining the feature that the physical force 
is still a quotient of quadratic expressions in this new parameter, as in the case of 
the original angular velocity parameter C and also for the physical component of the 
relative velocity in each point of view. Note that the physical force is a fractional 
quadratic function of the relative velocity, as it is also of any new variable differing 
from the relative velocity by a fractional linear transformation, like C or y, or the 
change in relative velocity from changing the observer. 

Comparing the full expression for the physical force expressed in terms of the 
relative acceleration, gravitoelectric, and gravitomagnetic terms in the threading 
point of view 

f{UY = 7(C/, m)2[K(0, mfviU, m)^^ - v{U, m)^H{mY - gimf] (4.5) 
leads to the identification 

g{mY = -K{(l>,mYv{U^,mfv{U+,mf, 

H{my = K{(f),mY[iy{U-,mf + iyiU+,mf] . (4.6) 

Similar relations hold for the hypersurface point of view with H{mY replaced by 
— 26'(n)r0 — H{n,e[)Y -\- 2k((/), n)'"j/(eo, n)'*. In the usual case in which the signed 
Lie curvature is negative, an upward gravitomagnetic vector field corresponds to a 
negative sum for the two geodesic velocities, as will be assumed in discussion below 
unless explicitly specified. 
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Consider the case of a circular orbit for which the Lie signed relative curvature 
K{4>,uy < 0, i.e., the familiar case in which the appropriate derivative of the Lie 
relative tangent vector is inward along the radial direction. Then the physical force 
is positive (outward) if i'{U- , u)'^ < 1/(11, u)'^ < y{U+,u)'^ and it is negative (inward) 
if v{U,u)^ > v{U+,u)* or v{U,uY < v{U-,uY. This result is exactly what one 
would expect in Newtonian gravity: namely, an outward push is necessary to remain 
on a circular orbit if the speed is less than the geodesic speed (Keplerian speed in 
the Newtonian context), while an inward push is necessary if the speed is larger 
than the geodesic speed. 

However, again as in Newtonian gravity, if the speed is larger than the geodesic 
speed and one increases it, the inward force necessary to maintain the circular orbit 
must also increase in magnitude and vice versa. This is not always true in the 
general relativistic case, not only in strong gravitational fields but also in certain 
weak gravitational fields. To investigate this property, one needs not only the sign 
of /{Uy but also the sign of its derivative with respect to the velocity of the test 
particle 1/(11, uy for a given point of view. 

This analysis is governed by the 3-parameter family of functions 

J^(iy; K, v-,v+) = k[v - - v+\/[l - v'^] (4.7) 

of the real variable v and by its derivative with respect to that variable. A specific 
spacetime and choice of observers collapses this family of functions to a 1-parameter 
family in which k,V-,v+ all depend on the radial coordinate. The way in which 
this collapse takes place leads to many possible behaviors. One must also be careful 
taking limits of this expression in which v and another velocity parameter go to ±1, 
since the order in which the limits are taken matters. 

The qualitative behavior of family of force versus velocity graphs (for different 
radii) is characterized first by the number of distinct real roots of the force function 
for a given radius, and second by the number of such roots which lie in the physical 
interval (—1,1) of "subluminal" values for nonzero rest mass test particle motion. 
One can distinguish three different regions of the spacetime (which may or may not 
be present) containing those radii for which one of the following conditions holds 
(assuming + vj^ < 0) 

region A: V- , both subluminal, 

region B: i/+ subluminal, z/_ superluminal, 

region C: both superluminal. 

When they exist, the interfaces between regions A and B and regions B and C are 
at the radii ?'(ab) ?'(bc) '^t which the circular null geodesies occur. 

Within each of these three regions, the qualitative properties of the force function 
are similar and follow from the expressions for its first and second derivatives 



dJ^(u; K, P-,p+)/di' 
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= -k[{u_ + U+){u' + 1) - 2v{l + - v'Y 

= + V+)[V^ - 2j//!/(,ei) + 1]/[1 - V^f 

= -k{v^ + V+)[v - Z/(crit)-]['^ - '^(crit)+]/[l - v"^? , 

d?'T{u', K, z/_, u^)/di''^ 

= -2k[{v^ + v+){v^ + - (3z/2 + 1)(1 + v-v+)]/[l - v^f 

= -2k{u. + + i)y - {iv^ + l)/i^(rei)]/[l - , (4.8) 



where 



i^(,ei) = {ly- + + V-V+) (4.9) 

up to sign is the relativistic difference of the two velocities and vanishes when 
i/_ + = as occurs in the Schwarzschild case. This is a subluminal velocity only 
when V- and z/+ are either both subluminal or both superluminal (regions A and 
C). 

The extrema of the force function are found by examining its critical points 
which only exist in regions A and C where i^(roi) is subluminal. The critical points 
occur at the roots (complex in region B) of the quadratic factor in its numerator 



Z^(crit)± = V±(l/(rel)) = 1 ± - V{j^\f /z^(rel) , (4-10) 

which satisfy 

J^(crit)-!^(crit)+ = 1 • (4-11) 

Under the assumption that < 0, they also satisfy z^(crit)+ < — 1 < J^(crit)- < 

in region A where —1 < i^(roi) < 0, and < z^(crit)- < 1 < i^(crit)+ in region C where 
< J^(rei) < 1- Thus the only critical point in the physical region occurs at the 
minus root 

l'(ext) = i^(crit)- = 7(rcl)i^(rcl)/(l + 7(rcl)) 

5'^(rci) {v^ + v^)l2 as i^(,.ci) . (4.12) 

In region A this leads to a maximum of the force function when k < (minimum 
when «; > 0) and a minimum in region C when «; < (maximum when k > 0), with 
the extreme value 

^(«^(ext);K, l'-,i^+) = + (1 + 7(rol))'^-i^+]/7(rol) • (4-13) 

In region B there are no critical points and no extrema, but one point of inflection. 
i/(].ei) has a vertical asymptote at V-V^^ = —1, dividing region B into two parts: on 
the region A side > —1) it is negative as in region A, while on the region C 

side {1^-1^+ < — 1) it is positive as in region C. 

In region A where both and 1^+ are subluminal, the value V(ext) has a simple 
and elegant interpretation in terms of 1-dimensional motion in special relativity. 
Consider three relative velocities , i^(cxt) of three observers with respect to a fixed 
fourth observer. Determine the velocity i'(ext) so that in its own rest frame, the 
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other two observers (corresponding to i'± ) have relative velocities which only differ 
in sign. Using the relativistic addition formula for velocities along a fixed direction, 
this condition is just 

- I^(oxt))/(l - l^+I^(0Xt)) = -{l^- - l'(0Xt))/(l - l^-i^CGXt)) ■ (4.14) 

Cross-multiplying this equation leads to the same quadratic equation in J^(cxt) as 



follows for u by setting the derivative (48) to zero. Thus the extremal force observer 



sees the two oppositely rotating geodesies with the same relative speed. 

It follows that the 4-velocity of the extremal velocity is just the renormalized 
average of the 4-velocities of the two oppositely rotating circular geodesies. This 
immediately yields the result 

;^(cxt) = I^(crit)- = [l-l"- + + 7+] ■ (4-15) 



From this and the reciprocal relation (4.11) one easily finds 



V{crit}+ = h-i'- - 7+i^+]/[7- - 7+] , (4-16) 

from which it follows that the 4-velocity of the second critical velocity is just the 
renormalized difference of the 4-velocities of the two oppositely rotating circular 
geodesies. 

The geometry of the relative observer plane of the circular motion (the t-tf) 
subspace of the tangent space) is very useful in visualizing the various velocities 
which arise here and in the spin precession analysis, and in extending this discussion 
to region C. For a given observer with 4-velocity u in this relative observer plane, 
denote the map which reflects across it by a tilde 

U ^ ^{u + lye) ^ U = ^{u - i^e) (4-17) 

and denote by a bar the commuting map which reflects (relative to U) across the 
nearest forward null direction to the orthogonal direction 

U = j{u + i^e) U ^ ji^{u+ 1/iye) , (4.18) 

where e is a unit vector orthogonal to u in the positive (f) direction in the relative 
observer plane of the circular motion. The pair {U, U) arises from either (u, e) or 
(m, — e) by a boost, for a timelike future-pointing U. For the threading observers. 



this corresponds to the bar map introduced in Eq. (2.12). Also let U± = 'y±{u+iy±e) 
be the geodesic (unit) 4-velocities whenever they are not null, whether timelike or 
spacelike, and the same for the critical 4-velocities [/(crit)± = 7(ciit)±(w -f J^(ciit)±e)- 
In regions A and C the latter pair are related to each other by the bar map because 
of Eq. dOTI ). 

Figure ^ illustrates the geometry of the various velocities in the relative observer 
plane of the circular motion for each of the three regions. In region A where the 
geodesic 4-velocities U± are both timelike, the relative difference velocity — i^(roi) > 
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is the (scalar) relative velocity of [/_ with respect to U+. The critical velocities cor- 
respond to the 4- velocities which are the normalized sum/difference of the geodesic 
4- velocities 

i^(crit)± = ^±{^{ic\}) , 

C/f„it)± = (C/^tC/:^)/||;7+tC/-|| . (4.19) 

In region C both U± are spacelike. Here one may rewrite the relative velocity 
^(rci) > in the following form 

I'M) = {i^--^ + + , (4.20) 

revealing that its sign-reversal may be interpreted as the relative velocity of tj - 
with respect to C/+, both timelike in region C. The critical velocities now correspond 
to the 4- velocities which are the normalized sum/ difference of the timelike barred 
geodesic 4-velocities or equivalently of the tachyonic four-velocities themselves 

c^rcrit)± = mTU^)i\\u+Ttj.\\ 

= {UI^UZ)/\\U+tU^\\ . (4.21) 

Thus in region C the extremal observer sees the two tachyonic geodesic velocities 
with equal but opposite velocities as in region A. 

In region B where J7+ is still timelike but [/_ is spacelike, the spin-critical 
velocities introduced below in section ^ are determined by a similar relationship 
but with z^(roi) replaced by its reciprocal (i.e., [/_ replaced by [/_) 

l'(crit,spin)± = V±(l/i^(i.cl)) = 1 ± y^l - l/i^(rol)^ '^(rol)- (4-22) 

But — t^(ioi) is the relative velocity of with respect to 11+ and — l/j^(ici) is the 
relative velocity of U- with respect to C/+. Therefore the corresponding spin-critical 
4-velocities are related to the critical ones in the same way that the critical ones 
are related to the geodesic ones 

f^(crit,spin)± = (^(CTit)+ T C^("rit)-)/ll^(crit)+ T t^(crit)-|| • (4.23) 

Thus if one reinterprets U± in Figure ||(6) as the critical velocities U(ciit) , then the 
four- velocities i7(sc) of the figure become the spin-critical velocities C/(crit,spin)- A 
similar diagram for the region B case < —1 has U+ and U- interchanged 

relative to u. 

In regions A and C one can introduce a family of "extremal force" observers with 
4-velocity f^^^t^ and its associated relative velocity V(oxt){u)'^ = i'(cxt)(t^(cxt), w)*^ for 
both the Kerr and Godel spacetimes, both of which have a pair of oppositely directed 



circular geodesies. The expression (4.13) with appropriate arguments then defines 



the acceleration of the extremal force observers. In region A their angular velocity 

C(cxt) = (r-0- + r+<^+)/(r_ + r+) < o (4.24) 
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Fig. 5. An illustration of the relative observer plane geometry of the various velocities which 
arise in the analysis of the physical radial force function for circular motion, for each of the three 
regions A, B, and C. In regions A and C, these determine the critical velocities (abbreviated by 
^(c)±) of the force function directly, while in region B where the critical velocities are complex, 
the geometry shown instead determines the spin-critical velocities (abbreviated by f7{s,c)±) for 
the spin precession function of section The case < — 1 of region B is obtained from the 

case —1<U-UJ^<0 shown in b) by interchanging U+ and U-, reversing the sign of the relative 
velocity. The case i/-u^ = —1 separating these two has instead U+ = U-. 



is easily evaluated by averaging the geodesic 4- velocities in the form (2.2), and is 
negative since r_ > r+ as explained at the beginning of section For low speeds 
(r± 1), this reduces to the average angular velocity C(gmp)- In region C their 
angular velocity is given by the analytic continuation of above formula (T± are 
purely imaginary so C(cxt) remains real) but is instead positive. 

Thus the extremal force condition picks out the family of observers which see 
two oppositely moving geodesic test particles symmetrically, i.e., with the same 
relative speed. Expressing the physical force function with respect to these observers 
leads to an expression which is an even function of v since ^(i^; k, — 1^+, = 
— — ly'^]. However, a pair of oppositely directed geodesic test particles 

which start from a given such observer do not return simultaneously to the same 
observer after each full revolution of the orbit because of the same asymmetry 



between corotation and counter-rotation which appears in the relationship (2.8) 
between the physical velocity component and the angular velocity. This is described 
below as a geodesic analog of the Sagnac effect. 
The difference 

C(ext) - C(gmp) = ^[(r- - r+)/(r_ + r+)][0_ - 0+] < o (4.25) 

is negative since r_ > r+ and < 0, 0+ > 0, so the extremal force observers 
counter-rotate with respect to the geodesic meeting point observers. The negative 
values of the angular velocities of both families of observers leading to counter- 
rotation with respect to the threading observers is due to the asymmetry in their 
equation of motion introduced by the gravitomagnetic vector force. In both cases 
these counter-rotate with respect to the slicing observers (which are determined by 
the null circular orbit meeting points) since the asymmetry introduced by the grav- 
itomagnetic force, being proportional to the speed, is smaller for timelike velocities 
compared to the speed of light, and so the corotating "dragging" effect is smaller. 

It seems reasonable that at least in the limit of small angular speeds, the physical 
force maintaining the orbits of the extremal force observers is maximum since in 
some sense they have the least relative angular motion with respect to the geometry, 
and hence require the most force to keep from falling towards the center of symmetry. 
This is intuitively clear in Kerr where the attraction towards the center is provided 
by the gravitoelectric field of the mass giving rise to the gravitational field. In 
the Godel spacetime in the slicing point of view, the gravitoelectric field is also 
inward, so a test particle initially at rest would also begin to fall towards the axis 
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of symmetry as it follows a circular orbit with a different center. 

Table 6. Signs of the physical component of the physical radial force and its velocity derivatives for 
circular orbits for the usual case K{<j>, uY < 0. The abbreviations v± = u{U±, and u = u{U, 
are used. 





Region A 


Region B 


Region C 




U- > -1 
u+ <l 


u- < -I 
v+ < 1 


ly- < -1 
iy+ > 1 


f{ur 


positive if : 
U- < v < u+ 
negative if : 
V > v+ ; 

U < V- 


positive if : 
u < u+ 
negative if : 

U > U-f- 


positive 


df(ur 

dv 


any sign 


ncgative 


any sign 




negative 


any sign 


positive 





Table | gives the sign of the physical force and of its first and second deriva- 
tives for each of these regions. For the Godel spacetime, there is always at least 
one timelike geodesic at every radius (the threading observers) so region C does 
not exist, while region C disappears from the Kerr spacetime in the extreme case 
a ~ A4. Similarly the region B, whose existence is due to the asymmetry in the 
angular motion introduced by the rotation of the reference frame, disappears in 
the Schwarzschild limit a = of Kerr. In the rotating Minkowski spacetime, the 
two real roots of the force function coincide, so there is only one circular geodesic 
(the counter-rotating slicing observer) and it is timelike at each radius where the 
threading point of view is valid, so only region A exists in both points of view. 
Figure || shows the dependence of the radius r as a function of a for Kerr for the 
two boundaries r'(AB), ''(bc) between the three regions A, B, and C, together with 
the radii of the relatively straight circles and embedding-signature-switching circles 
for the threading point of view discussed in the appendix. 

For each value of the radial variable, one can plot the physical force f{UY versus 
the physical component v{U, u)'^ of the relative velocity in the angular direction for 
both the threading and slicing points of view. Figure ^ shows a representative 
sample of these curves for the rotating Minkowski spacetime, the Godel spacetime, 
the Kerr spacetime with a/ Ai = 0.5, and the Schwarzschild spacetime a = 0. Shown 
also is the curve whose intersection with each force graph occurs at the velocity of 
the other family of observers. In the Godel slicing case, for example, this is just the 
positive horizontal axis since the theading observers are geodesic and corotating. 
The curves of relative extrema of the force graphs are also plotted. 

Note that in the slicing point of view in Minkowski spacetime, only the space 
curvature force due to the circular orbit is nonzero. This is just the usual inward 
centripetal acceleration required to maintain a circular orbit under the influence of 
no real forces, scaled by a squared gamma factor. It has a maximum at the relative 
velocity of the slicing observers, which are the extremal force observers in this case. 
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Increasing the radius flattens the force graph due to its inverse dependence on r. In 
the threading point of view, the force graph also migrates to the left as it follows 
the velocity of the single counter-rotating circular geodesic followed by the slicing 
observers at each radius, going infinite at the threading observer horizon radius. 

For a rotating spacetime the physical force f{Uy is not a symmetric function 
of the relative velocity y{U,u)'^ since the gravitomagnetic vector force introduces 
an asymmetry between the corotating and counter-rotating orbits. For the Kerr 
and Godel spacetimes the maxima and minima of the force graphs in regions A and 
C do not occur at zero velocity. At these cxtrcma the derivative of the force with 
respect to the velocity changes sign at a nonzero relative velocity, leading to effects 
which run counter to our ordinary Newtonian intuition about the radial force for 
circular orbits in a force field. 

For those spacetimes with a region A (Kerr and Godel), the sign of the radial 
force in that region behaves according to our Newtonian intuition about circular 
orbits. A positive (outward) radial force is required to support a circular orbit 
with speed less than the Keplerian speed, and a negative (inward) force is required 
to hold in a circular orbit with more than that speed. In the Schwarzschild case, 
the extrema of the force graphs occur at zero velocity and are local maxima in this 
region, so increasing the speed decreases the outward radial force needed to maintain 
the orbit. However, for the Kerr and Godel spacetimes, the local maxima occur at 
negative (slicing or threading) velocities, so for the counter-rotating orbits between 
zero velocity and the one at the local maximum, increasing the speed increases the 
outward radial force necessary to maintain the orbit, which is the opposite of what 
happens in Newtonian gravity. We refer to this simply as "counter-intuitive radial 
force behavior." Thus in region A, extending to infinite radius in Kerr but from zero 
radius out to the radius at which the counter-rotating geodesic becomes spacelike in 
the Godel spacetime, those counter-rotating orbits with speeds less than the speed 
at which the maximum of the force occurs experience this behavior. 

For those spacetimes with a region B (Kerr and Godel) where only the counter- 
rotating geodesic is spacelike (superluminal velocity), there is no local extremum 
on each force graph in this region. A particle on a counter-rotating circular orbit is 
forced to move with a velocity smaller than the (unphysical) geodesic speed and so it 
must be pushed outward (positive force) to maintain the orbit. This is as expected. 
However, if the speed of a counter-rotating orbit increases then the positive force 
increases as well, which is again counterintuitive. 

In the region C in the Kerr spacetime where both the corotating and the counter- 
rotating circular geodesies are spacelike (superluminal velocities), both corotating 
and counter-rotating test particles are forced to move with a speed less than the 
geodesic speed thus always yielding a positive (outward) radial force. As their speed 
approaches 1, the Lorentz gamma factor goes to infinity causing the physical force to 
increase to infinity as well. Since a local minimum occurs at a corotating (positive) 
velocity, increasing the speed increases the outward force for velocities outside the 
interval between zero velocity and the positive velocity where the minimum occurs. 
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Fig. 6. / versus u: plots of the radial physical component of the force acting on a nonzero rest 
mass test particle moving on an equatorial circular orbit as a function of the physical component 
of the relative velocity of the particle in the angular direction for the rotating Minkowski, Godel, 
Kerr {a/M = 0.5), and Schwarzschild (a = 0) spacetimes. The solid curves are the graphs of the 
force plotted versus the relative velocity for selected values of the radius. For the spacetimes other 
than rotating Minkowski, regions A, B, and C are the regions which are separated by the long 
dashed lines, from bottom to top. Region C is missing in Godel and region B in Schwarzschild. 
The short dashed lines connect the local extrema of the force graphs. The long-short dashed line 
marks the point on a given such graph which occurs at the velocity of the other family of observers. 
For Kerr the force curves move higher along the vertical axis with decreasing radius from region 
A extending out to infinity through region B into region C containing the ergosphere and finally 
approaching the observer horizon. This is also true for Schwarzschild except that region B dis- 
appears between the interface of regions A and C at r = 3m. For Godel they move higher along 
the vertical axis (slicing diagram) or the dashed-dotted line (threading diagram) with increasing 
radius from region A near r = into region B approaching the observer horizons. For the rotating 
Minkowski slicing diagram, they move outward from the vertical axis with increasing radius, and 
in the threading diagram they move to the left with increasing radius toward the observer horizon. 
The thickened force curve represents the radius of the ergosphere (threading observer horizon) in 
the Kerr slicing case and the corresponding slicing observer horizon radius in the Godel threading 
case. Almost no change in the force curve occurs beyond that radius in the Godel threading case. 



Thus some corotating and all counter-rotating orbits exhibit counter-intuitive radial 
force behavior. 

This happens for all circular orbits in region C in the Schwarzschild spacetime, 
where the interval between zero velocity and the velocity of this minimum shrinks 
to zero width. The corotating and counter-rotating geodesic velocities take the 
relativistically corrected Keplerian values 

v{U^,u)^ = TVM/r/y/l-2M/r (4.26) 

and the physical force reduces to 

f{Uf = K{cf>,uf^{U,ufMU,uf)^ - {,y{U^,um , (4.27) 

where the signed relative curvature is 

K{(f>, uY = - ^/l-2M/r/r . (4.28) 

The two circular geodesies become null at r — 3A4 where {^{Uzp, u^l = 1, and then 
spacelike for smaller radii. At this coordinate radius the physical force reduces to 
f{Uy = —K{(j), uY > which is independent of the velocity of the test particle, 
leading to the dashed line force graph separating the two regions A and C in Fig- 
ure ^(d). For 2A4 < r < 3A^ (region C), the circular geodesies are both spacelike 
and the physical force on the test particle must be positive. Moreover f{U) is an 
increasing function of |z^(J7, u)*^! in that region and so increasing the speed of the 
test particle increases the outward physical force, contrary to Newtonian intuition. 

Sis_^counter-intuitive radial force behavior was first discovered by Abramo- 
Ej in the Schwarzschild case inside r = 3M, the outer boundary in the 
equatorial plane of his "rotosphere" ,E3 and he explains it as due to a reversal of 
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the "optical centrifugal force" in a decomposition of the physical force into a sum 
of terms, de Felicec3'Ea has noticed this behavior in the Kerr spacetime using a 
factorization of the physical force in terms of angular velocity, which is closely 
connected to the threading point of view relative velocity, and he uses the term 
"prehorizon regime" to characterize it, since initially he interpreted this behavior 
as a way of signaling the approach to a horizon. Later he discovered that it can 
occur for counter-rotating orbits even a great distances from a rotating objected 
Barrabes, Boisseau, and Israelahave used the hypersurface point of view to describe 
the same phenomenon for both Schwarzschild and Kerr, but using a factorization 
of the physical force in terms of the slicing relative velocities. 

Of course all of these effects depend on the zero point used for the relative speed, 
and are therefore observer-dependent. The choice of observer used to describe the 
effect depends on exactly what one wants to measure or explain. The threading 
point of view is relevant to behavior as seen from infinity, while the hypersurface 
point of view might be appropriate for local considerations like accretion disks. For 
the extremal force observers in region A introduced above, the extrema of the phys- 
ical force occur at zero relative velocity. A similar family exists in region C. For 
example, in the region C of the Kerr diagram, one sees that the corotating family 
of observers for which the minimum occurs approaches the slicing observers at the 
horizon but increasingly corotates faster than those observers until it experiences 
its own outer observer horizon as one moves out to the radius of the corotating 
photon orbit. In this region C the counter- intuitive radial force behavior occurs for 
all relative velocities with respect to these observers just as in the Schwarzschild 
case. For a certain interval of larger radii (region B), no extremum occurs. Then 
in the region A outside the radius at which the counter-rotating photon orbit oc- 
curs, one has a counter-rotating family for which the maxima occur. This second 
family has an inner horizon at that radius and counter-rotates with respect to the 
threading observers but approaches them at infinite radius. For this second family 
the counter-intuitive radial force behavior does not occur. Figure 0(a) shows the 
coordinate angular velocity of each of the geometrically defined observer families for 
the Kerr snacetime as a function of the radius (see Fig. 2 of de Felice and Usseglio- 
Tomassetc3), including the Carter family associated with the usual orthogonal frame 
in which the Kerr metric is^ated in Boyer-Lindquist coordinates (Eq. (33.2) of Mis- 
ner, Thorne, and WheeleiS), an observer family which has a coordinate angular 
velocity C(car) = a/(r^ + a^) and simultaneously diagonalizes the electric and mag- 
netic parts of the curvature tensor. Figure ^(6) shows the corresponding plot of the 
observers' physical relative velocities in the slicing point of view. 

Finally consider the less usual case in which the Lie signed relative curvature 
is positive, K{(f>,uY > 0, which would seem to reverse the sign of the physical 
force as one crosses the Lie relatively straight trajectory from the side for which 
k{4>,uY < 0. However, the superluminal velocity i'{U-,u) blows up at this radius, 
switching signs so that the physical force remains finite and does not change sign, 
even though the relative centripetal acceleration itself does change sign ( "reversal of 
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the centrifugal force" in the language of Abramowicz et al,Ei but in the true relative 
geometry rather than the optical one). For example, in the extreme Kerr case in 
the threading point of view, the relative curvature becomes positive for r < r(i.s) 
and the space curvature force along the corotating circular orbit becomes negative 
while v{U-^ mY blows up there. In the Godel spacetime this instead occurs in the 
slicing point of view where the relative curvature becomes positive for r > r(i.s) • 

5. Optical geometry and inertial forces in the static case 

The preceding companion article [BCJl] shows how the optical gauge conformal 
transformation of the spatial metric changes the relative centripetal acceleration to 
the optical relative centripetal acceleration and reshuffles the various spatial forces 
in the spatial equation of motion. As discussed there, the optical geometry is only 
natural in the static case, which in the present context applies to the Schwarzschild 
limit of the Kerr spacetime. 

Those results are easily specialized to this latter case, so that one can see how the 
various radial spatial forces for circular orbits are transformed under the conformal 
transformation. Agreeing that orthonormal components of optical quantities are 
normalized in the optical geometry, the transformation of the signed Lie relative 
curvature (equation (14.3) of [BCJl]) is 

k{(t),uY = cr-Y'i{(t),uY - (Incr),?] , (5.1) 

while the curvatures themselves are 

k{(l),uY = -{l-3M/r)/r , k{(I>,uY ^ -{I ~ 2M/rY^^ /r . (5.2) 

With the static relation g{uY = (Incr) ^ that follows from the optical choice a = 
or a = respectively, this leads to equation (A. 23) of [BCJl] for the trans- 
formation of the centripetal acceleration, here specialized to orthonormal compo- 
nents 

= <^tUu,uY HU,uY\'g{uY . (5.3) 

For comparison with the literature, note that the conformal factor is not present in 
the covariant (unnormalized) coordinate component form of this equation. 



Using this result in equation (4.2), one finds the following expression for the 
physical force necessary to maintain the circular orbit in terms of the Lie optical 
relative acceleration 

fiUY = <J-'liU, uYa\^^^^{U, uY - l{U, u)fI;^]{U, uY + l{U, ufHU, uf\''g{uY 
= a-'^{U,uf~af^^{U,uY-g{uY, (5.4) 

where the extra term 7(f/, M)^|i^(f/, u)''^|^g(u)'' which comes from the conformal 
transformation of the centripetal acceleration combines with the term — j{U, uYg{uY 
in the spatial gravitational force to form a velocity-independent force, namely mi- 
nus the gravitoelectric field —g{uY- Abramowicz refers to this as the gravitational 
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Fig. 7. (a). (J versus f: plots of the coordinate angular velocities of the preferred observer families 
in the equatorial plane of the Kerr spacetime as a function of the radial coordinate f for a = 0.5. 
The thick outer long dashed curves above, below) correspond to the null orbits which 
enclose the physical region of timelike orbits, a region divided in two by the thick short dashed 
curve C(sli) ; these three curves meet at the horizon at the left edge of the graph. The thick long- 
short dashed horizontal axis corresponds to C(thd) = 0- The two solid curves immediately above 
and below C(sli) ^^'^ C{thd) = respectively are C(car) ^i^d C{gmp)- The thin dashed curves are 
C(goo)± = 0±i which intersect the null orbit curves at the interface between regions A, B, and C 
by definition. The solid curves in regions A and C terminating in these same intersection points 
correspond to (C{ext))- The curve decreasing from the upper geodesic photon point to the lower 
geodesic photon point in region B is C{crit,spin) where the threading precession is extremal. 

(b) . V versus f: the plot of the corresponding physical relative velocities in the slicing point of 
view as a function of the radius. 

(c) . V versus f: the same plot as in (6) showing the three physical relative velocities at which the 
corrected spin precession (^{U,m, E)" vanishes. 



forcejHy but with respect to the local rest space of the test particle itself, not the 
static observers fixed in the geometry of the spacetime. 

The remaining term in the physical force, namely the sign reversal of the Lie 
optical relative centripetal acceleration multiplied by the proper time correction 
factor 7(C/, m)^ and the conformal transformation factor for relative acceleration, 
is what Abramowicz identifies as his centrifugal force in the static case. The Lie 
optical relative centripetal acceleration is (see [BCJl]) is explicitly 



"(lie) 



iU,u) 



|i>(C/,j.)'^|2[lng(«)^^ 



1/21 



(5.5) 



In their attcmr 



|S to extend Abramowicz's work to 
uate this result for the slicing 



where v{U,uY = (jv{U^u) 

the Kerr spacetime, Iyer and Prasannall3'C::acn e 
decomposition, while Prasanna and ChakrabartillS evaluate it for the threading 
decomposition, modulo conformal factors. 

The step which goes from the spatial equation of motion expressed in terms of the 
sum of the relative centripetal acceleration and the spatial gravitoelectric force to 
the same equation expressed in terms of the sum of the optical relative centripetal 
acceleration and the Abramowicz gravitational force in the Schwarzschild case is 
explicitly 



f{Uf 



-i{U,u)MU,uf? 2 



r-iM 



(1 - 2M/r)-i/2 
(l~2M/r)~^'^ . (5.6) 



This clearly shows the result that at r = the physical force acting on the test 
particle is independent of its velocity. This last property of the physical force is 
lost in the corresponding stationary case because of the presence of the velocity- 
dependent gravitomagnetic force. 
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6. Spin precession 

A test gyroscope following a test particle world line in spacetime undergoes 
Fermi- Walker transport along that worldline. Precession of the spin direction is 
a relative effect, depending on both the family of test observers and the single 
test observer following the test particle world line. Spin precession for geodesic 
motion will be analyzed first, and then followed by the case of accelerated motion. 
PapapetrouB has shown that a spinning particle does not follow a geodesic in 
spacetime but is subject to Riemannian tensor forces. These forces and any non- 
gravitational forces acting on a gyroscope are represented by the force F{U,u)°', 
which leads to the Thomas precession term given below for accelerated motion. 

The precession angular velocity of the spin of a gyroscope (with 4- velocity U") 
as seen in its own local rest space is a relative effect, parametrizing the derivative 
of the relative boost between the local rest space of an observer congruence and the 
local rest space of the gyroscope along its world line. Given an observer-adapted 
orthonormal spatial frame {Ea"} on spacetime which is tied to the observer con- 
gruence by spatial co-rotating Fermi- Walker transport, one can boost the frame to 
the local rest space of the gyroscope, and the components of its spin with respect 
to that boosted frame will undergo a time-dependent rotation with an associated 
proper time angular velocity having components 

C(C/,u,i?)''=7(C/,u)[C(cfw)(C/,")" + C(sc)(C/,",£^)"] , (6.1) 

where 

C(cfw)(c/,w)" = -^H{ur 

-j{U, u)[7(f/, u) + l]-i[z^((7, u) x„ F{U, u)r 

= C(gm)(^)"+C(thom)(C/,wr+C(gco)(f/,")" (6.2) 

is the precession angular velocity relative to a spatial (with respect to u") co- 
rotating Fermi- Walker frame transported along U°' (and consisting of gravitomag- 
netic, Thomas, and geodesic terms), while the space curvature term 

C(.c)(C/,u,i?)'^ = ^v{ur'^Tiu)[t\d\c]HU,uf (6.3) 

is the relative angular velocity of the latter frame with respect to the observer- 
adapted frame {Ea"}, and the factor "/{U, u) in Eq. ( |6.1| ) is necessary to convert the 
observer proper time derivative to the gyroscope proper time derivative. Conversely, 



removing the gamma factor in Eq. (6.1) one obtains the precession of the spin as seen 
by the sequence of observers along its path, removed of the directional distortions 
caused by the Lorentz boost between the two local rest spaces. 

This description is valid for an arbitrary spacetime.El For the threading observers 
in the nonlinear reference frame associated with post-Newtonian coordinates in that 
approximation to general relativity, one obtains the Schiff formula for the precession 
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with respect to the "distant stars." One can use those same general results in the 
threading point of view for any of the three spacetimes under consideration to 
obtain an exact precession formula for a gyroscope in an equatorial plane circular 
orbit. The natural observer-adapted (right-handed) orthonormal spatial frame 
which is tied to the threading observer congruence by corotating Fermi- Walker 
transport (in the equatorial plane) consists of unit vectors along the coordinate 
radial r direction, the observer local rest space </> direction, and the —6 or positive z 
direction, respectively. The relative observer boost of this spatial frame to the local 
rest space of the gyro, completed to a spacetime frame by the fivro 4-velocity, is 
called a phase-locked frame by de Felice and Usseglio-TomassetU'Ei who investigate 
who investigate the Kerr case, and it is also the spacetime Serret-Frenet frame for 
the circular orbit discussed in detail by Iyer and VishveshwaraJ^^ who treat the 
general stationary axisymmetric case. 

For Kerr, the resulting spin precession formula has the same interpretation as 
the Schiff formula, namely the precession of the spin with respect to a local frame 
locked onto the distant stars. For the Godel spacetime, one obtains the precession 
with respect to a local frame locked into the perfect fluid source of the gravitational 
field, which reduces to the gravitomagnetic term alone in the corotating geodesic 
case where the relative velocity is zero. For the rotating Minkowski spacetime, 
the circular geodesies correspond to points fixed in the inertial coordinates of the 
slicing observers, so the precession question is not interesting without considering 
accelerated circular orbits, where the Thomas precession is obtained in the slicing 
point of viewB In all cases the only nonvanishing component of the spin precession 
is 

mm.Ef = 7(f/,m)[C(efw)(C/,m)^ + C(sc)(t/,"i,i?)'] 

= 7(f/,"i)C(«ph) ((7, m,i?)^ (6.4) 

where C(sph)(C^, i5)^ is the precession of the spin with respect to the static spher- 
ical frame as seen by the sequence of observers along the gyro world line. 

First consider the case of geodesic motion so that the Thomas precession is 
zero. For an equatorial circular geodesic in an axisymmetric stationary spacetime, 
the angular velocity terms in the threading point of view take the form 



(6.5) 



where 
and 



Fl^\{U,mY = ^{U,mMmY + ij.([/, m)^i/(m)^] , (6.6) 

C(sc)iU,m,Ey - -ii/(;7,m)^(ln7^^),, = i.(C/,m)'^«(0,m)^ 

= iy{U,mY sgn{K{(j),mY) R{m)/p{(/),m) , (6.7) 



where R{m) = 
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For the Kerr spacetime the relative velocities and gamma factors of the circular 
geodesies are 

i^{U±,m)^ = VA/[a±^/r/M{r-2M)] , (6.8) 

2M. 



(1 



-)(1 



-1/2 



and the gravitomagnetic, geodesic, and space curvature angular velocities are 

C(gm)(TO) = 
C(geo)(C^±,m)'' = 



= ± 



C{sc){U±,m,Ey = - 



MjiU±, m)L'{U±, m)^{-VA + aiyiU±,m)^) 

{j{U±,m) + iy{r - 2M) 
7(j7±,m) - 1 
liU±,m) V ' 

p{U±, m)^[r{r ~ 2Mf - Ma"^] 



aM 



r^{r-2M) 



Ar^{r-2M) 



(6.9) 



The sum of these three terms 



C{sph){U±,m,Ey = C(sm){mY +C{eeo){U±,mf + C^^^){U±,m,Ey 

= Tl{U±,m)-'^ (6.10) 

gives the total precession angular velocity of the spin with respect to the spherical 
static frame in terms of the threading observer proper time. It vanishes in the limit 
'y{U±,m)~^ — > as \iy{U±,m)\ — > 1 which occurs for the two freefall photon orbits, 
corresponding to a locking of the spin to that frame. This is consistent with the 
picture of a photon as a massless spinning particle with spin along the direction of 
motion, i.e., locked to the 0-direction in the case of a circular orbit. On the other 
hand, the precession from the point of view of the gyro has the same simple form 
C{U±,m, Ey = ^yTM/r^ as in the Schwarzschild Hmit. 

However, since the static spherical frame is locked to the radial direction, it 
undergoes a rotation of 27rsgnr^(J7, n)'^ during each revolution with respect to 
Cartesian-like frames along the gyro worldline and so the spin rotates an addi- 
tional — 27rsgni^(J7, n)"^ with respect to the Cartesian-like frames. This corresponds 
to a precession angular velocity of the spherical frame components of the spin of 
the gyroscope equal to the orbital coordinate angular velocity which must therefore 
be added to the space curvature precession to correct it, onc;c suitably scaled to 
account for the difference between coordinate and observer proper times 



d(p±/dT(^u±,m) = '^{U±,m)'^ = ± 



1 



2M 



1 



2M 



±a 



(6.11) 
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The corrected space curvature angular velocity is then 

C(sc,cor) 



= [1 + sgn{K{(f>, mf)R{m)l p{U±,m)]v{U±,mY 



(6.12) 



aM 



r^{r-2M) 




1 - 



2M 



1 



2M 



±a^ 



For r > r(i.s) when sgn(«;((^, m)'') < 0, this reduces to [A(m)/27r]i^(J7±, m)"^, where 
the quantity the quantity A(m) is the deficit angle of the tangent cone to the embed- 
ding of the r-(j) coordinate surface (with the threading geometry) in E3 , as explained 
in the appendix. Thus Thome's conical deficit argument using the threading em- 
bedding spaceO holds exactly in this region. 

The total precession angular velocity of the gyroscope with respect to the distant 
stars can then be written 



7(C/±, m)[C(gm) {my + C(gco) {U±, my- + C(sc,cor)(C^±, m, Ey] 
T^/?[l-(l-^±2a,/5)-/T (6.13) 



Figure^ shows the gravitomagnetic, geodesic, and corrected space curvature angular 
velocity for extreme Kerr as functions of r/A4. 

Table 7. The spin precession <^(f7, m)^ terms for geodesic circular orbits, including low speed limits 
indicated by arrows. The abbreviation 7_ = ■y{U—,m) is used. Note that 7_ = C(l — S^)~^^^ in 
the Godel case. 





Rotating 
Minkowski : U— 


Godel: 
U- 


Kerr : 

u± 


''(gm) 




-n 


_^ _ aM 


^(gco) 





-n[i - 7i^] ^ 




'{sc,cor) 


^7- [7_ - 1] ^ 


2n[i - i/c] -» 




^(cor) 


-Q 


n[i - 2/Vi - 52] ^ -n 


r 1 3 a/M , / 



Table summarizes the corresponding results for the rotating Minkowski and 
Godel spacetimes, including their small r (low relative velocity) limits. This ta- 
ble also lists the small Ai, small a limits for Kerr for comparison. These latter 
expressions agree with the well known post-Newtonian results, but it is impor- 
tant to note that the geodesic term in the precession formula and the corrected 
space curvature term, arising from very different phenomena in the strong field 
case, combine to form the single "geodetic precession" term in this weak field limit. 
For the counter-rotating circular geodesies in rotating Minkowski spacetime, the 
precession corresponds to the rotation of the inertial coordinate axes with respect 
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(a) Kerr: counter- rotating geo 
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(b) Kerr: corotating geo 
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Fig. 8. C versus f: plots of the geodesic, space curvature (corrected), and gravitomagnetic preces- 
sion angular velocity of the spin of a gyroscope moving on (a) counter-rotating and (b) corotating 
equatorial circular geodesies of a Kerr black hole with a = 0.5 as functions of the radius f . Their 
sum gives the corrected precession Qj-q^j ((7, m, i?)^/7([/, m) with respect to the threading ob- 
servers, shown as the long dashed curve. Since these are all finite, the corrected precession terms 
themselves all diverge at the ergosphere boundary. 
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to the threading observer axes. For the corotating circular geodesies in the Godel 
spacetime, only the gravitomagnetic precession is nonzero and has the same value 
as for the counter- rotating geodesies. 

Th^a^results agree with those of Rindler and PerlickEl and Iyer and Vishvesh- 
wara,t§E3 both of whom use rotating coordinates t' = t, (j)' = (f) — ujt for which the 
circular orbit has a fixed angular coordinate in order to find the precession angle 
of the gyroscope after one loop. This precession angle Acj) and the precession 
angular velocity C(cor)(^^±j EY are related by a scale factor and an integral over 
the loop 



= ±[2V^^(^7±,m)^]7(^7±,m)-lC(cor)(C/±,m,^;)^ (6.14) 

using equation ( |6.1l|) and recalling the proper time conversion factor. Comparing 
this with equation ( p. 12 ) and recalling the additional gamma factor in equation 
(5.1) shows that the contribution to the total precession angle from the corrected 
space curvature precession equals the threading deficit angle A(m) modulo the sign 

A(/.±(sc,cor) =±A(m) . (6.15) 

Of course, with some more formula juggling one can obtain the spin precession 
for the case of circular orbits with arbitrary acceleration. To express this as a 
function of the relative velocity and radius, one need only replace F{U, m,Y in the 
radial component of the u = m version of Eq. ( |6.2| ) by ^{U, m)~'^ f{U, mY , where the 
physical force f{U,mY is itself expressed in terms of those variables by Eq. ( |4.3| ). 
For example, the Thomas precession term has the form 

Citi.ora){U,mY - u{U,mY[j{U,m) + l]-'f{U,mY 

= K{(j), mYHU, mY[-f{U, to) + l]"^7(t/, mY 

X [iy{U, mY ~ i^{U-,mY] [iy{U, mY - v{U+, m)^] . (6.16) 

Gathering the preliminary expressions for the other precession terms and doing 
some involved algebra in which Eqs. (4.6) are helpful leads to the following general 
result for the total precession angular velocity 

C{U,m,EY = ^-/{U,m)-^{dT/d:^){iy{U,mY;K{(j),mY,iy{U^,mY,iy{U+,mY) , 
C(cor)(C/,m,^)^ =C(C/,m,^)^+7(C/,m)z.(;7,TO)V7^//2 , (6.17) 

generalizing Eq. ( |6.10 ) to the accelerated case for all the spacetimes under consid- 
eration. 

The precession angular velocity ({U, m, BY thus vanishes at the critical points 
of the radial physical force function, corresponding to the locking of the spin to 
the static spherical frame for the extremal force observers. These observers act as 
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the boundary between the counter-rotating and corotating spin precession orbits 
relative to the radial direction. Table || correlates the sign of this precession with 
these intervals of relative velocity values in each of the three regions A, B, and C 
assuming < 0. In the limiting Schwarzschild case where region B disappears 

and !^(ext) — 0; the spin precession ({U,m,Ey has the opposite sign compared to 
the angular velocity in region A and the same sign in region C. 

Table 8. Signs of the spin precession for circular orbits of a given angular relative velocity. 





Region A 


Region B 


Region C 


C(U,m,E)^ > 


(-l,!^(oxt)) 




(i^(cxt), 1) 


C(U,m,E)^ = 






!^(cxt) 


C(U,m,E)^ < 


('^(cxt), 1) 


(-1,1) 


(-l.'^Ccxt)) 



Using the properties of the function !F introduced in Eq. (4.7) and the relation 



(4.6), the gyro spin precession angular velocity with respect to the static spherical 
frame from its own point of view takes a very simple form analogous to the physical 
acceleration it experiences 

a{UY = T{v{U,mY]K{4),mY,v{U-,mf,v{U+,mf). (6.18) 

One can then carry over to the spin precession the entire analysis of the extrema of 
the physical force function a(UY . Defining a relative difference velocity by Eq. (4.9) 
with iy± replaced by i/(crit)± leads to i^(roi,spin) = l/i^(ici) which is in the physical 
range only in region B where |i^(rci) I > 1; ^(crit)± ^-re complex, and the spin precession 
is always negative. 

Only in this region are the corresponding critical values related to this new 
relative velocity by Eq. (4.10) real 



^(crit,spin)± = V±(l/t^(rcl)) = i^(rcl) 1 ± 1 ^ ^ / '^{ic\) 



(6.19) 



As before only the minus root is in the physical range leading to the spin extremal 
velocity 

' ' (6.20) 



I', 



(cxt,spin) 



^(rcl) 



(rol) 



at which a maximum negative (but minimum absolute value) spin precession angular 
velocity occurs with value 

- ^{'^{cxt,spm)',H/2, I^(crit)-: 1^(001) + ) = -^^[2 + 7(rcKspin)] /7(rcl,spin) • (6.21) 

Figure H(6) illustrates the geometry of the relative observer plane associated with 
this analysis. 
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Figure || shows both (a) the spin precession C,{U,m,EY and (6) the corrected 
spin precession C(coi ) {U, m, EY at the same selected radii as a function of the thread- 
ing relative velocity for a Kerr black hole with a — 0.5 outside the ergosphere where 
the spin precession with respect to the distant stars makes sense. The intersections 
of these curves with the short dashed curves in Figure |9|(6) occur at the velocities 
of the geodesic orbits, corresponding to Figure ||. The long dashed curves mark 
the boundaries r(AB)j ^(bc) between the regions A, B, and C. One can see three 
families of observers in region A for which the corrected precession vanishes. Upon 
closer inspection one sees that there is one pair which respectively corotate and 
counterrotate faster than the pair of timelike circular geodesies and share the same 
respective observer horizons, and another family near the slicing observers which 
has the ergosphere as its horizon. 

Working instead with the coordinate angular velocity, de Felicelil has shown 
(following from a comparison of Refs. |3l| and (28|, see also Ref. |2j) and implied 
the fact that the corotating Fermi- Walker spin precession vanishes exactly at the 
two critical points of the radial force function, and is posilive between them and 
negative outside. In the approach of Iyer and VishveshwaraES for stationary axially 
symmetric spacetimes, specialized to their equatorial plane case T2 — Q — A[2) , the 
precession angular velocity ({U,m,E)^ is just the sign-reversal of the first Serret- 
Frenet torsion ti and the acceleration a{UY of the test particle (physical force) is 
the Serret-Frenet curvature k. For the equatorial plane case their Eq. (50) in their 
notation reduces to 

n = ^A/^3d^K = ~l[MTiU,mYj^/^^Y'dca{UY 

= -^j{U,m)-^d,J^ = -C{U,m,EY . (6.22) 

This first torsion is called the Fermi drag by de Felice. In the case of general motion 
along a Killing trajectory in a stationary axially symmetric spacetime, one can show 
that the gyro spin precession angular velocity (Eq. (20) of Ref. |3^) relative to the 
Frenet-Serret frame E'q" is just 

- c.(;7)(FS)" - P{U, uYpCiU, u, eY ■ (6.23) 

7. Sagnac effect and synchronization defect 

The Sagnac effectiiii and its timelike analog, both of which in turn are con- 
nected to the synchronization defect ,c3 refer to the asymmetry in the arrival times 
of a pair of oppositely rotating timelike geodesic or null circular orbits at a given 
radius as seen by a given rotating observer. If (Ci,C2) is the ordered pair of coor- 
dinate angular velocities of such a pair (either (C_,C+) or (</)_, 04.)), and C is the 
angular velocity of a rotating observer with 4- velocity U distinct from this pair, one 
easily finds that the difference and average of the coordinate arrival times after one 
complete revolution with respect to this observer are 



At = 5(C;Ci,C2) = i2-ti = 2^[l/(C2-C)-l/(C-Ci)] 
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Fig. 9. f versus u: plots of (a) the precession angular velocity C,(U,m, EY and (b) the corrected 
precession angular velocity C{cor)('^i "'■i °f spin of a gyroscope versus the orbital physical 
angular velocity for selected radii for equatorial circular orbits of a Kerr black hole with a = 0.5. 
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= -47r[C - (Ci + C2)/2]/[(C - Ci)(C - C2)] , 

i(avg) = (tl+t2)/2-2^[l/(C2-C) + l/(C-Cl)] 

= ^[Ci-C2]/[(C-Ci)(C-C2)] , (7.1) 
the ratio of which gives the relative difference in the coordinate arrival times 

AiA(avg) = 4[C - (Ci + C2)/2]/[(Ci - C2)/2] . (7.2) 
For the pair of oppositely rotating timelike geodesies one has 

At(geo)(t/) = 5(C;^_,0+) = -4^[C-C(gmp)]/[(C-0-)(C-0+)] , (7.3) 
while for the pair of oppositely rotating null orbits one has 

At(„un)(C/)=5(C;C-,C+) = -47r[C-C(nmp)]/[(C-C-)(C-C+)] , (7.4) 

recalling that C(nmp) — C(si)- Each of these may be evaluated for the slicing, thread- 
ing, and extremal force observers, when the three orbits are distinct for a given 
application. Note that the null arrival time difference is proportional to the angular 
momentum (2.11) of U. 

The Sagnac time difference and its timelike geodesic analog are given when the 
observer is taken to be the threading observer 

At(nuU)(m) = 5(0;C-,C+) =4K(nmp)/(C-C+) 

= 47r(C--i+C+"')/2 = 47rM^ , 
At(geo)(m) = 5(O;</«_,0+) = 47rC(gmp)/(0-'/>+) 

= 47r(0_-i +0+-i)/2 = 47ra [Kerr]. (7.5) 

The null difference is positive/negative when the threading observers corotate/- 
counter-rotate with respect to the slicing observers (sgn A^jnun) = sgn N"^ = sgn M^). 
The geodesic difference is positive assuming an upward gravitomagnetic field 
(sgnAt(gco) — sgniJ^), with the counterrotating geodesic returning first to the 
threading observer {t2 < ti, where 2 labels the corotating geodesic). The geodesic 
formula is not valid for Godel since the threading observer is itself geodesic, nor for 
Minkowski which only has one geodesic. 

The synchronization defect is just half the Sagnac time difference and indicates 
the change in coordinate time which occurs during one corotating spatial loop (</> 
coordinate line in spacetime) 

/■27r 

Ai(SD)(m)=/ M^d(j) = 2TTM4, = At(^nuii){m)/2 . (7.6) 

The analogous coordinate time difference for the geodesic case in the Kerr spacetime 

27r/C(car) = Ai(goo)(TO)/2 = 27ra (7.7) 
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corresponds exactly to one threading loop of a circular curve which is spatial with 
respect to a Carter observer (C(car) = a/(a^ + r^), so C(car) = 1/a by Eq. ( 2.12| )), 



connecting the average-time-of-return point on the threading observer world line 
with either return point on the same world line. 

Figure ^ shows the relationship between the arrival times of oppositely rotating 
circular geodesic and null orbits and the Sagnac effect and its timelike geodesic 
analog in region A. This is illustrated for the case in which the shift is negative as in 
the Kerr spacetime, i.e., the threading observers counter-rotate with respect to the 
slicing observers, so the corotating null orbit returns first to the threading observer 
from which it originated. Independent of the sign of the shift but assuming that 
the gravitomagnetic field is upward, the geodesic meeting point observer counter- 
rotates with respect to both the slicing and threading observers, and the extremal 
force observers counter-rotate in turn with respect to the geodesic meeting point 
observers. Even though the extremal force observers see the oppositely rotating 
pair of geodesies moving with the same speed, the counter-rotating geodesic arrives 
after the corotating one as shown. 

The horizontal unit vector e{n) in Figure ^ is along the slicing relative veloc- 
ity direction tangent to the coordinate line (constant t). The horizontal short 
dashed line connecting P± contains all the (horizontal) slicing relative velocities. 
The unit vector e(m) is along the threading relative velocity direction and indicates 
the threading local rest space angular direction. Extending this forward around 
one loop of the cylinder leads to the change in the coordinate time equal to the 
synchronization defect. 

One can express the threading observer synchronization defect, Sagnac effect, 
and timelike geodesic generalization of it using a single formula involving the phys- 
ical fields. The 1-form 

-f{U,m)-^Uo,dx" = [ma, + l/i^{U,mfe{m)o,]dx" 

= -M{dt - M^dcj)) + 700i/Vzy(C/, mf d4) (7.8) 

orthogonal to the 4- velocity C/" restricts to zero along the world line, so the observer 



Fig. 10. The conformal diagram of the front half [— 7r/2, it/2] x R of the flattened i-fli coordinate 
cylinder for a fixed r in region A, showing the null and geodesic time differences for corotating and 
counter-rotating orbits for the various observer families. The vertical axis is along the (null meeting 
point) slicing observers (n), the time coordinate axis is along the threading observer direction (m), 
and the horizontal axis is along the ip coordinate direction (constant t). The diagram illustrates 
the case in which the slicing observers corotate with respect to the threading observers as in Kerr 
so that the time lines tilt to the left. The relative observer plane tangent space (region A case) is 
superimposed on the cylinder together with the corresponding world lines emanating from a given 
single initial point. The slicing, threading, geodesic meeting point, and extremal force observer 
4-velocities and world lines are shown, together with the null 4-momenta = it e{n)°' of 
the two oppositely rotating null orbits. The synchronization defect is shown at the bottom of the 
figure corresponding to one loop of a threading purely spatial curve with unit tangent e(m)". 
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proper time arrival time is 

T{m) = Mt{m)^[sgniy{U,m)'f'] [MM^ + -/^^'''^ /v{U,mf] 







- 27rAf/|C| . (7.9) 
Thus difference in arrival times of two oppositely directed world lines is then 

Ar(m) = M[t2{m) - tx{m)] = 2n \mM^ + -i^^^f^[l/v{Ui,mf + l/v{U2,mf] . 

(7.10) 

For the oppositely directed null orbits {v{P, m)"^ = il)i this just gives 

AT(„un)(™) = 47rMA/^ = 47rp^(m) , (7.11) 

namely the conserved angular momentum of the threading observers. For the op- 
positely directed timelike geodesic orbits this just gives 

AT(gco) (m) = 47rp^(m) - C{m)H{mY / g{my , (7.12) 

1 /2 

where C{m) = 27r7^^ is the threading circumference of a (/> coordinate circle, using 
Eq. (|^) to re-express the quotient of the sum and product of the geodesic relative 
velocities. For the limit v{Ui,m) — > ±cxd in which the trajectories become spatial 
curves with respect to m, this gives twice the synchronization defect (defined to be 
the change in time around one such loop, not two) which has the same value as the 
proper time Sagnac effect 

AT(gD)(TO) = 5AT(„un)(m) . (7.13) 

These same considerations should apply to any stationary family of observers, 
once they are expressed in terms of the threading potentials and fields associated 
with the family. For the slicing observers, for example, both the angular momentum 
and the threading gravitomagnetic field vanish so the Sagnac and synchronization 
defects vanish. 

Similarly the slicing point of view expression for the proper time geodesic effect 
has an analogous formula with the substitutions H{mY — > ~29{n)^., C{m) 
C'in) = 2^5^^i/2 

Ar(geo) (n) = + 2C{72)9{n)^'/giny , (7.14) 
since the slicing angular momentum vanishes. 

8. Concluding Remarks 

Test particle behavior is one of the most important tools in interpreting the local 
spacetime structure of a given spacetime, whether it is a nice exact solution with 
symmetry which may be worked with analytically or an approximate solution with 
no symmetry derived with very complicated mathematical machinery. The three 
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familiar spacetimes analyzed here in terms of natural test observer families each 
have a rich and interesting geometry that weaves together in slightly different ways 
the various effects that their gravitational fields have on the motion and spin of test 
particles. 

By adopting a general approach to the test observer splitting and the geometry 
they use to interpret their measurements, many diverse issues heretofore described 
in widely differing languages are brought together in a single discussion based on 
familiar observer-defined quantities. Instead of vague references to centrifugal and 
Coriolis forces that motivate the interpretation of fields in linearized approximations 
to general relativity or the somewhat obscure formal representations of such forces 
in the fully nonlinear theory, this general approach gives a clear and precise rep- 
resentation of the many ways these concepts can manifest themselves in the arena 
of Einstein's theory. By analyzing these three familiar spacetimes in this approach, 
and showing the power that it has in interpreting physical aspects of their geome- 
try and showing the wider applicability of the analysis to allow comparisons among 
them, we hope to encourage others to consider its use when appropriate. 
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A. Embedding cross-section diagrams 

For the 2-dimensional metrics in polar-like coordinates {r, (j)} under considera- 
tion 

hrrdr'^ + h^4,d(l}'^ = hnndB? + R^dcj)^ , (A.l) 

where the components depend only on the radial coordinate, the embedding into 

or with cylindrical coordinates {R, 4>, Z} can be accomplished in two steps. 
First one inverts the transformation from the radial coordinate r to the radial 
circumferential arclength coordinate R = h^^'^/'^, for which the circumference of 
a radial coordinate circle has the usual expression 2'kR. Then one performs an 
integral to get the curve in the R-Z half plane which is the constant (j) cross-section 
of the embedding surface. It is convenient for graphing purposes to introduce the 
rescaled variables R = R/TZ and Z = Z/TZ analogous to f = r/TZ, where TZ is the 
length scale of Table II. The case Hrr = 1 corresponds to the flat case of E2 in polar 
coordinates. When him > one next identifies the result with the flat cylindrical 
coordinate metric on E3 restricted to the surface Z = Z{R) 

dR^ + R'^dcj)^ + dZ{Rf = [l + {dZ{R) /dRf]dR^ + R^dcj)^ . (A.2) 

If instead hun < 1, one identifies the result with the flat cylindrical coordinate 
metric on M3 restricted to the surface Z = Z{R), where Z is now a timelike 
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coordinate 

dR^ + R^dcf)^ - dZ{Rf = [1 - {dZ{R)/dRf]dR^ + R^dcf)"^ . (A.3) 

In both cases this gives a first order differential equation for Z[R) with an integral 
solution as in the Schwarzschild caseEa 

1 ± {dZ{R)/dRf = IiRR = hrr{dr/dRf , (A.4) 

or equivalently 

\dZ{R)/dR\ = ^±{hRB. - 1) , 



\Z{R)-Zo\ - 




(A.5) 



where the plus/minus sign corresponds to the embedding in E-^/M^. Dropping the 
absolute value sign, Z may be chosen initially to increase with i?, but a negative 
derivative may later be required if the embedding turns back towards the Z-axis as 
Z increases. 

In general it is not possible to express the argument of the square root explicitly 
as a function of R but one can make a parametric plot of Z{r) and R{r) using 

\dZ{r)/dr\ - ^±[hrr-{dR{r)/drY] , 

R{r) = VK^, (A.6) 

which has the integral solution 

\Z{r) -Zo\= f ^±[hrr-{dR{r)/drf]dr . (A.7) 

This may be numerically integrated if necessary. 

For cases with a regular origin i? — > as r — > 0, the initial conditions Zq = and 
ro — puts the origin of the {r, cj)} coordinate system at the origin of the embedding 
space. In other cases, one may start at a physically interesting value of r for rg, 
which then determines R, and choose Zq conveniently. 

Once the curve in the R-Z plane is obtained either analytically or numerically, 
the full surface is obtained by revolving it around the Z axis. Introducing the Lie 
relative radius of curvature from [BCJl] for this context 

p{4,,u) = \k{4,,uY\-' = I - (lni?),,//i„i/2|-i (A.8) 

1/2 

leads to the result /i^'^ — p{(f>, u)/R and 

\dZ{R)/dR\ ^ v/±[p(0,w)Vi?2 - 1] = v/±[p(0,M)2-i?2]/i? . (A.9) 

The ratio p{(j),u)/R is identically 1 for a flat space, greater than 1 for a Euclidean 
embedding, and less than 1 for a Minkowski embedding. In the latter cases, if 
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this ratio passes through 1 at some radius, the tangent to the R-Z cross-section is 
horizontal and the embedding space signature switches sign. This tangent hne is 
instead vertical when \dZ[R)/dR\ —> oo, which occurs when p((j),u) — > cxd at a Lie 
relatively straight circle. On either side of the radius r(ss) at which the signature 
switches sign, the choice of sign for the derivative dZ{r)/dr is free. This sign will 
be always be chosen to be positive to join the two embedding spaces on either side 
of a half plane, in which case 

sgn(K((/), uY) = - sgn{dZ/dR) . (A. 10) 

Embedding diagrams for which the signature changes have been used by Smarr0 
in studying the 2-surface cross-section of the horizon of a charged Kerr black hole. 

When the tangent to the embedding cross-section is not vertical, one can intro- 
duce the tangent cone obtained by revolving the tangent line at a point {R, Z{R)) 



on the cross-section. The relation (A. 9) then shows that p{(j),u) must be identified 
with the distance of the tangent line point of contact from the vertex of the cone 
calculated in the embedding space geometry, thus giving a nice geometrical inter- 
pretation to the Lie relative radius of curvature. For the £'3 case, the defect angle 
of the cone is easily found to have the positive value 

A = 27r[l-i?M0,u)] , (A.ll) 

which in the usual case in which sgn(K((/), to)'') < determines the space curva- 
ture contribution to the precession angle of the spin of a gyroscope follow ing the 
circular orbit after one revolution modulo the sign, as discussed by Thorne.Ell The 
same formula holds for the spacelike cones in M3, where the defect angle is instead 
negative. Note that in approaching a relatively straight trajectory (vertical tangent 
in the R-Z cross-section) this defect angle approaches 27r, while in approaching a 
horizontal tangent in the R-Z cross-section the angle approaches the flat value 0. 

The Gaussian curvature of this surface is easily evaluated from standard formulas 
to be 

K = ~{Rhrr^/'^)-^{R,r/hrr^/'^),r = -R^^R^ff . (A.12) 

For the rotating Minkowski slicing and the Godel threading cases where hrr — 1, 
this reduces to K = —R~^R rr- This may also be rewritten as 

K = sgn{K{(j),uY) R-^[R/p{(l},u)]^f. . (A.13) 

Thus the sign of the Gaussian curvature depends on the signs of the radial relative 
centripetal acceleration (sign of —R^r) and of the radial derivative of the ratio 
R/p{(j),u). It is also straightforward to show that 

sgn{K) = Tsgnid^R/dZ^) , (A.14) 



which implies that when the cross-section curve is concave away/toward the Z-axis 
in E3, the Gaussian curvature is negative/positive, and vice versa in M3. 
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A.l. Rotating Minkowski spacetime 

The spatial geometry in the shcing point of view is fiat, so the usual Euclidean 
geometry applies and the Lie signed relative curvature of the coordinate lines is 
just the reciprocal of the radial coordinate r. 

The r-(f) surface in the threading point of view is instead an inhomogeneous 
surface of revolution of negative Gaussian curvature —2>^^/Ti? which starts at the 
value —2>/Ti? at the origin and as r increases it decreases to — oo at the observer 
horizon at the value f(h) = 1- The embedding is in M3 and the singularity at the 
observer horizon occurs dX R —> 00. Figure ^l](a) shows the constant cj) cross-section 
of this surface, with the horizontal and vertical axes showing the rescaled variables 
TZ and Z respectively. 

The threading case 

Since dR/dr = 7'^ is always positive, then hnn = < 1 and the embedding 
occurs in M3 with a regular origin and a single sign for dZ{R)/dR which may be 
taken to be positive. The relationship R = (with 7 = (1 — f2)-i/2) inverts 
to f = 7~i^ (with 7 = (1 + R'^y/'^). Dropping the absolute value signs in the 



embedding equation (A. 5) and choosing Zq = leads to 

Z{R) ^ / - {l + x^)-^dx . (A.15) 
The change of variable w ~ 2\n^ = ln{l + R^) = — ln(l — f^) leads to 



i / ^yl + 2coshydy , (A.16) 
Jo 



which in turn suggests the further change of variable u — sinh(w/2) = (7 — 7 ^)/2 
and 

Z{u) = 2 /" V(3/4 + i'2)/(l + w2)dv 
"'0 

= -|[F(a, 1/2) - F{Tr/2, 1/2)] + 2[£;(«, 1/2) - E{7r/2, 1/2)] 



+2mv/(1 + "')/(3/4 + w2) , (A.17) 

where a = arctan(l/M) hjL,a slightly corrected version of formula 3.169.2 on p. 276 of 
Gradshteyn and Ryzhik,o and F{(p, k) and E{ip, k) are elliptic integrals of the first 
and second kind in their notation. (When using MAPLE V and Mathematica re- 
spectively, the correspondence is i?(</3, k) — LegendreE(sin93, k) = EllipticE((p, fc^).) 



A. 2. Godel spacetime 

Figure ^ shows the constant </> cross-section of the r-ip surface with the hori- 
zontal and vertical axes showing the rescaled variables TZ and Z respectively in the 
(b) threading and (c) slicing points of view. 



Intrinsic Derivatives and Centrifugal Forces in General Relativity 



53 



(a) Godel: Sli 




1 



Fig. 11. Z versus R: embedding cross-section diagram for the r-<p surfaces in the threading point 
of view for the Godel (a) slicing and (fe) threading points of view and (c) the rotating Minkowski 
threading point of view. In the latter two cases the asymptotic null lines in the M3 embedding are 
shown with respective intercepts Z — 0.567f (iior) ^^^d ^ ~ ~0-406f (iior) which are comparable 
fractions of the horizon scale, but which appear much different since units of vorticity rather than 
curvature are used in the graph. In the first case starting at the origin with a Euclidean embedding, 
the embedding surface moves away from the vertical axis and then turns back at the radius of the 
single relatively straight circle, next switching signature and finally asymptotically approaching 
the vertical axis again at the observer horizon along an asymptotic null cone. 
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In the threading point of view the spatial metric of the r-(j) surface is in one of 
the standard forms for a constant negative-curvature geometry {K = — \/7i?) and is 
easily embedded in A'l^i as a spacelike pseudosphere of radius TZ. The circumference 
R of the circles of constant radial coordinate increases without bound but the mag- 
nitude of the Lie signed relative curvature approaches a lower limit corresponding 
to a maximum Lie relative radius of curvature of magnitude 7?., a consequence of 
the negative curvature geometry. 

In the slicing point of view where r is a radial arclength coordinate this surface 
has positive Gaussian curvature K = (27?.^)~^(l-t-2s^— 2s'*) / (1— s^)^ which increases 
from the value 1/ {2Ti?) at the origin to infinity at the observer horizon. R begins to 
decrease with increasing r after reaching the Lie relatively straight circle at the value 
r(rs) = 2arcsinh(2~*/^) w 1.317 (-R(rs) — 1) where the Lie signed relative curvature 
changes sign and the embedding turns back toward the Z-axis. Beyond that the Lie 
relative centripetal acceleration points outward, i.e., towards larger r values (but 
smaller R values) as occurs for the optical relative centripetal acceleration close 



= V2V2 - 2 K, 0.910, 0.847). Finally the embedding surface 

returns to the Z-axis (i? — 0) at the value r^^) — 2arcsinh(l) w 1.763 [Z = IZ^y^) 
corresponding to the observer horizon, where the embedding is singular. 

The threading case 

Here the standard embedding with r = at the origin of A/3 is a pseudosphere 
of radius TZ 



The slicing case 

The origin r = is again regular, so Zq = is appropriate, and one may 
drop the absolute value signs initially in the explicit embedding equation ( [A.5| ). 
Since dR/dr — c(l — 2s^)/-\/l — changes sign as r increases, this embedding 
turns back toward the Z-axis, after which dR/dr passes through the value 1 where 
the embedding switches from to A4^. Attempting to invert the relationship 
R{u) = 2-\/u(l — u) = 2s\n~^s^ by solving a quadratic equation for u ~ s'^ 



in terms of which (r — 2 arcsinh(\/M)). However, i? is a doublevalued function of r 
or u, so the latter variables are more appropriate to use. The integral 




Z = -Tl+ \JR^ + TV- . 



(A.18) 




(A.19) 
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to 



2y/{l/2-v^)/{l-v^)dv (A.20) 
can be integrated in terms of elliptic integrals depending on the range of u. 
Embedding in E^. Starting at the regular origin with u e [0, -^] this evaluates 

Z = 2E (^a,l/V2^ - F (^a,l/V2^ , (A.21) 

where a = arcsin(-\/2u), by formula (3.169.9) on p. 276 of Gradshteyn and RyzhikH 
Note that the endpoint value u = l/\/2 at the signature change corresponds to 
a — tt/2, while the relatively straight circle aX u = 1/2 corresponds to a = 7r/4. 
Let = 2E{tt/2, 1/V2) - F{n /2,1/V2) ~ 0.847. 

• Embedding in M3. The new initial condition at m = 1/V2 for u G [l/\/2, 1] 
with the change in sign of the square root argument leads to 



Z - = 2E{(3, 1/V2) - F{I3, 1/V2) - 2^/{u^ -l/2){l-u^)/u , (A.22) 

where B = arcsin (V2u2 - 1/u), by formula (3.169.11) on p. 277 of Gradshteyn and 
Ryzhik.L3 The observer horizon then occurs as the limiting point where u — 1, 
R^O and P = n/2, so that Z = 2%^). 

A. 3. The Kerr spacetime 

The Kerr case can only be handled analytically in the limit of the Schwarzschild 
spacetime. The embedding geometry of the r-cf) equatorial plane in the Schwarzschild 
spacetime is well known to be a parabola of revolution in E^ with its vertex at the 
observer horizon which occurs at the relatively straight circle ( "throat" ) .H The 



integral (A. 7) leads to 

Z{f) = y^2/{u-2)du = 2V2{f - 2)^/2 , (A.23) 



with Zq = and fo = 2 while R = f, explicitly confirming the parabolic embedding 
cross-section. 

The Kerr equatorial plane must be studied numerically in both points of view 
since inverting the relationship between r and the circumferential coordinate R or 



performing the integral (A. 7) cannot be done analytically. The slicing embedding 
diagrams have been given by Sharptj while a useful discussion of the general problem 
has been given by Romano and Price.LJ 

For the slicing point of view the plus sign is always valid and the embedding 
occurs in E^. The initial conditions are taken to he Zq = and tq = ^(h) = 
1 + %/l - e [1, 2], leading to ^(r(h)) = 2. Note that R"^ =f^+a'^ + 2c?- jf implies 
dR/dr = {f^ — a?)/{r'^R) > 0, which is positive and finite for the physical ranges 
< a < 1 and r > r{h)- Since A = at r = r(h), then hrr — > 00 as r — > f'(h)^7 so 



56 D. Bini, P. Carini and R. T. Jantzen 



that also dZ/dR = [dZ{r)/dr\ / [dR{r) / dr] oo. Thus all the slicing embedding 
cross-sections start out at ^ = 2, Z = in the R-Z plane with a vertical tangent as 
in the Schwarzschild case where that point is the vertex of a parabola. 



For the threading point of view the integrand of (A. 7) switches sign, while as 
f — > ^(erg)^ — 2+, _R — > oo, so the observer horizon is pushed out to infinity in M3 
as in the Minkowski threading case. Therefore initial conditions are taken to be 
Zq = and fo is chosen to be the single real zero of the integrand which occurs for 
the physical ranges 0<a<l,f>2, equivalent to 

2f2 {f-2f + d^ [2f{f - 2)2 - a^] = . (A.24) 

As shown in Figure these roots lie in the interval [2,2.344] (approximately) with 
the endpoint values respectively corresponding to a = and a = 1. The example 
d = 0.5 has 7~(ss) w 2.157. 

For r > tq the plus sign is relevant and the embedding is in , while for f < fg 
the minus sign is relevant and the embedding is in M3. Thus the i?-axis is taken 
to be the signature change line for all values of a, with above and M3 below. 
One also finds that dR{r)/dr = at the single real solution of f{f — 2)^ ~ d^ — 
for the allowed values or f and a, leading to a vertical tangent in the R-Z plane 
at the Lie relatively straight circle, where the embedding cross-section turns away 
from the Z-axis. 

The Gaussian curvature for both the threading and slicing cases is negative 

n [f(f - 2)2 + 3a2(f- 1)1 ^ 
■A^'^(th) = 4(r--2)^ ' ^^-2^) 

(-1) - f4[f(f2 + a2)-h2a2]2 ' 

with the common Schwarzschild limit A4^K = —1/f^. For Kerr it becomes infinite 
at the threading observer horizon as in the rotating Minkowski case, while remaining 
finite at the slicing observer horizon (where it vanishes for d — 1, representing 
a relatively straight circle). Figure [ij shows the constant cross-section of the 
embedding surface for different values of the rotation parameter a in both points of 
view. 
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Corrections 

This reformatted version contains the following misprint corrections of the orig- 
inal article (to which the page numbers refer) and one reference publication update: 

• p. 146. Table 1 Line 1: Godel lapse N change Vl — to Vl — 

• p. 147. Eq. (2.3) Line 1: change second minus sign to an equal sign 

• p. 150 Table 2. Row for signed relative curvature, Line 2, Godel (last) column: 
equal sign missing after kappa symbol before minus sign 

• p. 169. Line following Eq. (4.25): change </> to ^ in both occurrences 

• p. 174. Figure 7.(a),(b): The steep downward curve in region B was replaced 
and the last line of figure caption (a) changed from: "The steep nearly straight 
curve in region B is ..." 

to 

"The curve decreasing from the upper geodesic photon point to the lower 
geodesic photon point in region B is ". . ." 

C(crit,spin) whcrc the threading precession is extremal. 

• p. 176. Eq. (5.4): second subscript "lie" should be in roman type 

• p. 182. Eq. (6.16) Line 1: remove 7^ factor 

• p. 188. Eq. (7.10): change 47r to 27r 

• p. 189. Eq. (7.12): change 2C(m) to C(m) 



